REPORT  DOCUMENTATION  PAGE 


AFRL-SR-AR-TR-07. 


Public  reporting  burden  for  Ihis  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  se 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  the 
this  burden  to  Department  of  Defense.  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports  (0704-0188).  1215  J 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
valid  OMB  control  number  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS 


0078 


1.  REPORT  DATE  (DD-MM-YYYY) 
30-01-2007 


2.  REPORT  TYPE 

Final  Technical 


3.  DATES  COVEKtu  irium  . 
01-11-2003  -  31-10-2006 


4.  TITLE  AND  SUBTITLE 

(U)  Laser  Diagnostics  for  Reacting  Flows 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

FA9550-  0*  -1-0009 


5c.  PROGRAM  ELEMENT  NUMBER 

61102F 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 

2308 


Ronald  K.  Hanson 


5e.  TASK  NUMBER 

BX 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESSES) 


Stanford  University 
651  Serra  Street 
Stanford  CA  94305 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESSEES) 

AFOSR/NA 

875  North  Randolph  Street 

Suite  325,  Room  3112 

Arlington  VA  22203-1768  tr,i  - 

J>r  rru\iajr\  -n 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Advanced  optical  diagnostic  techniques  relevant  to  propulsion  were  investigated.  The  techniques  studied  were  based  on  laser 
spectroscopy,  with  emphasis  on  spectrally-resolved  absorption  and  laser-induced  fluorescence  (LIF).  Laser  sources  included  tunable  cw 
near-infrared  diode  lasers  and  tunable  (or  fixed- wavelength)  pulsed  lasers  operated  at  ultraviolet  (UV)  or  infrared  (I R)  wavelengths.  The 
cw  lasers  were  spectrally  narrow,  allowing  study  of  innovative  diagnostics  based  on  spectral  lineshapes,  while  the  pulsed  lasers  provided 
intense  bursts  of  photons  needed  for  techniques  based  on  LIF.  Accomplishments  of  note  included:  (1)  development  of  a  new  imaging 
diagnostic  based  on  infrared  planar  laser-induced  fluorescence  (IR  PLIF),  (2)  investigations  of  quantitative  ultraviolet  (UV)  PLIF  of  NO 
and  C02  in  high-pressure  combustion  environments,  (3)  the  development  of  a  new  temperature  diagnostic  using  UV  absorption  of  C02  for 
high-temperature  combustion  environments,  (4)  development  of  advanced  wavelength-multiplexed  diode  laser  absorption  sensing  of  non- 
uniform  temperature  distributions,  gas  temperature  in  scramjet  flows,  and  tunable  mid-IR-based  fuel  sensing,  and  (5)  further  development 
of  quantitative  tracers  to  image  fuel  distribution  using  ketones  and  the  aromatic  toluene.  The  full  spectrum  of  results  was  published  in 
thirty-eight  papers  in  the  AIAA  and  peer  reviewed  literature,  seven  PhD  theses,  and  forty-three  presentations  and  invited  lectures. 


15.  SUBJECT  TERMS 

Laser  Diagnostics,  Combustion,  Propulsion,  Absorption,  Fluorescence,  Temperature,  High  Pressure 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

Unclassified 


b.  ABSTRACT 

Unclassified 


c.  THIS  PAGE 

Unclassified 


17.  LIMITATION 
OF  ABSTRACT 


UL 


18.  NUMBER 
OF  PAGES 


72 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Julian  M.  Tishkoff 


19b.  TELEPHONE  NUMBER  (indude  area 
code) 

(703)696-8478 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


^STANFORD 

•  mechanical 


engineering 


Final  Technical  Report 


LASER  DIAGNOSTICS  FOR  REACTING 
FLOWS 


Grant  AFOSR  FA9550-04-0009 


distribution  statement  a 

Approved  for  Public  Release 
Distribution  Unlimited 


Prepared  for 

Air  Force  Office  of  Scientific  Research 
Dr.  Julian  Tishkoff,  Technical  Monitor 


For  the  Period 

November  1,  2003  to  October  31,  2006 


Submitted  by 

Professor  Ronald  K.  Hanson,  Principal  Investigator 


HIGH  TEMPERATURE  GASDYNAMICS  LABORATORY 
Mechanical  Engineering  Department 
Stanford  University,  Stanford  CA  94305 


Table  of  Contents 


1.0  Introduction . 4 

2.0  Project  Summaries . 5 

2.1.  Infrared  PLIF  imaging  diagnostics . 6 

2.2  Quantitative  LIF  and  PLIF  measurements  in  high  pressure-flames . 13 

2.2.1  Strategies  for  quantitative  NO  LIF . 13 

2.2.2  First  PLIF  imaging  of  CO2 . 19 

2.3  Temperature  diagnostics  based  on  CO2  absorption  of  UV  light . 24 

2.4  Advanced  wavelength-multiplexed  diode  laser  absorption  diagnostics . 30 

2.4.1  Non-uniform  temperature  distributions . 30 

2.4.2  Scanned-wavelength  modulation  for  high  SNR  absorption . 42 

2.4.3  Mid-IR  measurements  of  fuels  using  diode  lasers . 52 

2.5  Photophysics  studies  in  support  of  tracers  for  fuel  imaging . 57 

3.0  Publications  and  Presentations . 61 

4.0  Personnel . 68 

5.0  Technology  Transitions/Transfers . 69 

5.1.  Interactions  and  collaboration  with  AFRL  . 69 

5.2.  Transitions . 70 

8.0  Inventions  or  Patent  Disclosures . 72 

9.0  Honors/Awards . 72 


3 


1.0  Introduction 

Advanced  optical  diagnostic  techniques  relevant  to  propulsion  are  investigated. 
The  techniques  studied  are  based  on  laser  spectroscopy,  with  emphasis  on  spectrally- 
resolved  absorption  and  laser-induced  fluorescence  (LIF).  Laser  sources  include  tunable 
cw  near-infrared  diode  lasers  and  tunable  (or  fixed-wavelength)  pulsed  lasers  operated  at 
ultraviolet  (UV)  or  infrared  (IR)  wavelengths.  The  cw  lasers  are  spectrally  narrow, 
allowing  study  of  innovative  diagnostics  based  on  spectral  lineshapes,  while  the  pulsed 
lasers  provide  intense  bursts  of  photons  needed  for  techniques  based  on  LIF. 
Accomplishments  of  note  include:  (1)  development  of  a  new  imaging  diagnostic  based 
on  infrared  planar  laser-induced  fluorescence  (IR  PLIF),  (2)  investigations  of  quantitative 
ultraviolet  (UV)  PLIF  of  NO  and  CO2  in  high-pressure  combustion  environments,  (3)  the 
development  of  a  new  temperature  diagnostic  using  UV  absorption  of  CO2  for  high- 
temperature  combustion  environments,  (4)  development  of  advanced  wavelength- 
multiplexed  diode  laser  absorption  sensing  of  non-uniform  temperature  distributions,  gas 
temperature  in  scramjet  flows,  and  tunable  mid-IR-based  fuel  sensing,  and  (5)  further 
development  of  quantitative  tracers  to  image  fuel  distribution  using  ketones  and  the 
aromatic  toluene.  The  full  spectrum  of  results  is  published  in  thirty-eight  papers  in  the 
AIAA  and  peer  reviewed  literature,  seven  PhD  theses,  and  forty-three  presentations  and 
invited  lectures. 
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2.0  Project  Summaries 


2.1  Infrared  PLIF  imaging  diagnostics 

Infrared  planar  laser-induced  fluorescence  (IRPLIF)  of  carbon  dioxide  (CO2)  has 
been  demonstrated  in  the  past  [2. 1.1 -4;  e.g.  see  Refs.  1-4  at  the  end  of  section  2.1]  for 
imaging  room  temperature  and  elevated  temperature  jets.  The  technique  has  also  been 
used  for  imaging  the  fuel  region  of  a  laminar  CO  diffusion  flame,  and  the  nascent 
formation  of  CO2  in  a  laminar  methane  diffusion  flame  [2. 1.2,3].  In  the  current  project, 
we  extend  this  technique  for  the  first  time  to  supersonic  flows,  using  IRPLIF  of  CO2  to 
image  the  exit  region  of  an  underexpanded  jet. 

Background 

IRPLIF  involves  the  excitation  of  molecular  vibrational  modes  via  infrared  laser 
radiation  (an  excitation  wavelength  of  2.0  pm  is  used  for  this  work),  followed  by  imaging 
of  the  subsequent  vibrational  fluorescence  (fluorescence  is  collected  from  the  4.3  pm 
band  of  C02).  Continued  improvement  in  infrared  laser  sources  and  infrared  camera 
technology  make  IRPLIF  a  competing  method  with  UV/visible  tracer  PLIF  for 
visualization  of  flows.  IRPLIF  offers  the  advantage,  in  some  instances,  of  not  needing  to 
add  a  tracer  to  the  flow  if  molecules  with  an  IR-active  vibrational  mode,  such  as  CO2,  are 
naturally  present  in  sufficient  abundance.  Even  in  situations  where  they  need  to  be  added, 
small  molecules  with  IR-active  vibrational  modes  offer  the  advantage  of  being  gaseous 
under  room  temperature  conditions.  This  reduces  experimental  complexity  by  doing 
away  with  the  need  to  seed  vapors  from  a  liquid  into  a  gaseous  flow.  Additionally,  for 
applications  in  air-based  systems,  these  molecules  have  small  impact  on  the  properties  of 
the  fluid  being  visualized  because  their  molecular  weights  and  specific  heats  are  close  to 
those  of  molecular  oxygen,  molecular  nitrogen,  and  water. 

The  sub-atmospheric  pressure  regime,  which  is  particularly  important  for  high-speed 
flows,  offers  additional  advantages  for  IRPLIF.  Electronic  PLIF  can  be  dominated  by 
collisionless  nonradiative  decay  processes  [2.1.5],  however,  for  IRPLIF  of  small 
molecules,  e.g.  CO2,  the  only  nonradiative  deactivation  processes  present  are  due  to 
collisions.  Because  the  deactivation  processes  involve  collisions,  the  deactivation  rates 
are  linearly  proportional  to  the  number  density  of  the  colliding  species.  Reducing  the 
pressure  of  the  environment  to  be  imaged  reduces  the  deactivation  rate  of  the 
fluorescence,  and  increases  the  fluorescence  quantum  yield.  For  high-speed  flows  this 
advantage  may  be  limited  by  the  tolerable  exposure  time  for  the  flow  of  interest.  Typical 
deactivation  times  for  vibrational  fluorescence  are  on  the  order  of  microseconds  for 
ambient  temperature  and  pressure  conditions.  By  contrast,  lifetimes  of  nanoseconds  or 
picoseconds  are  typical  for  electronic  fluorescence  so  the  motion  of  the  flow  is  easily 
frozen.  IRPLIF  imaging  in  supersonic  and  hypersonic  applications  may  thus  require  the 
sacrifice  of  signal  level  to  achieve  short  integration  times  necessary  to  freeze  the  flow 
motion.  The  following  sub-sections  present  IRPLIF  theory  and  analysis  as  well  as  the 
experimental  setup  and  imaging  results  for  an  underexpanded  jet  of  CO2. 
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1RPLIF  Theory  and  Analysis 

The  fluorescence  equation  in  the  linear  excitation  limit  can  be  written  as[2.1.5] 


(2.1.1) 


where  Sf  is  the  number  of  photons  per  pixel  incident  on  the  camera,  E  is  the  laser  fluence 
[J/cm2],  h  is  Planck’s  constant  [J-s],  v0  is  the  center  frequency  of  the  laser  [1/s],  dVc  is  the 
volume  of  the  pixel  in  the  laser  sheet  [cm3],  xabs  is  the  mole  fraction  of  the  species  being 
excited,  P  is  the  local  pressure  [MPa],  k  is  the  Boltzmann  constant  [J/K],  T  is  the  local 
temperature  [K],  g(v)  is  the  laser  lineshape  function  [s],  a(P,T,  g(v))  is  the  absorption 
cross  section  of  the  species  being  excited  [cm2],  <f>  is  the  fluorescence  quantum  yield 
(FQY)  which  depends  upon  the  local  pressure,  temperature,  and  species  mole  fractions 
Cf»)»  and  rjco/i  is  the  optical  collection  efficiency  which  includes  transmission  efficiencies 
of  the  optics  and  transmission  of  the  fluorescence  through  the  gas  (radiative  trapping 
influence)  as  well  as  the  collection  efficiency  of  the  optics. 

Quantities  of  key  interest  in  Eqn.  (2.1.1)  are  the  absorbing  species  mole  fraction,  the 
local  pressure,  and  the  local  temperature.  The  absorption  cross  section  and  FQY  are 
important  because  they  are  dependent  upon  the  quantities  of  interest  and  influence  the 
signal  level  of  the  measurements.  In  supersonic  flows,  the  pressure,  temperature,  and 
mole  fraction  are  nonuniform,  necessitating  modeling  of  the  functional  dependence  of  the 
absorption  cross  section  and  FQY  to  determine  the  dependence  of  the  fluorescence  on 
these  quantities. 

Model  Results 

The  absorption  cross  section  is  dependent  upon  temperature  through  the  Boltzmann 
fraction  in  the  rotational  line  being  excited,  as  well  as  through  the  temperature 
dependence  of  the  absorption  lineshape  due  to  collision  broadening  and  Doppler 
broadening.  The  pressure  dependence  of  the  cross  section  arises  due  to  the  influence  of 
collision  broadening.  In  addition,  to  these  influences,  the  cross  section  depends  upon  the 
convolution  of  the  absorption  and  laser  lineshapes.  Absorption  cross  sections  are 
calculated  using  line  parameters  from  the  HITRAN  database  and  fitting  them  with  Voigt 
profiles.  These  lineshapes  are  then  convolved  with  the  laser  lineshape  which  is  modeled 
as  Lorentzian  with  a  FWHM  of  500  MHz  based  on  the  manufacturer’s  specifications. 
Cross  sections  are  calculated  over  a  two-dimensional  range  of  temperature  and  pressure 
for  multiple  possible  excitation  wavelengths.  The  plots  of  absorption  cross  section  are 
then  used  as  an  aid  in  selecting  excitation  wavelengths. 

The  definition  of  FQY  for  IRPLIF  differs  slightly  from  that  for  electronic  PLIF  due  to 
the  longer  fluorescence  lifetimes  associated  with  infrared  fluorescence.  For  IRPLIF  the 
fluorescence  lifetime  is  often  longer  than  the  camera  integration  time  so  the  FQY  is 
defined  as  an  integration  of  the  spontaneous  emission  rate,  from  the  states  emitting  into 
the  wavelength  collection  band  of  the  camera,  normalized  by  the  total  number  of  photons 
absorbed. 

The  deactivation  of  the  excited  vibrational  population  is  due  almost  completely  to 
collisions,  although  a  small  contribution  is  due  to  the  fluorescence  itself.  For  the  small 
molecules  that  we  are  treating  here,  collisions  are  the  mechanism  by  which  vibrational 
energy  is  transferred  between  molecules  (intermolecular)  and  between  modes  of  a 
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molecule  (intramolecular).  Assuming  translational  and  rotational  equilibrium  the 
problem  reduces  to  finding  the  time-dependent  population  of  each  vibrational  state  for  the 
molecules  in  the  mixture  of  interest.  The  approach  taken  utilizes  the  master  equation 
formalism  where  the  population  evolution  in  each  state  has  the  following  form  [2.1.6] 

=  K Wjntra  +  .inter  +  ^W.sk  +  ^ VT  +  +  ^  pump  (2-  1  -2) 

where  nv  s  is  the  time-dependent  population  in  vibrational  state  v  of  molecule  s,  Kw.imra  is 
the  sum  of  intramodal  vibration-to-vibration  (VV)  transfer  rates  into  the  vibrational  state 
v,  Kyy jmer  is  the  sum  of  the  intermoda!  VV  transfer  rates  into  the  vibrational  state  v,  Kyy,,k 
is  the  sum  of  the  intermolecular  VV  transfer  rates  into  the  vibrational  state  v  for  collision 
between  molecules  of  type  s  and  molecules  of  type  k ,  Kyj  is  the  sum  of  the  vibration-to- 
translation  (VT)  relaxation  rates,  KAu  is  the  sum  of  the  spontaneous  emission  rates  into 
state  v,  and  Kpump  is  the  laser  pumping  rate  into  the  state  v,  which  is  zero  for  states  not 
coupled  by  the  laser  radiation.  To  solve  for  the  time  evolution  of  population  in  each  of 
the  vibrational  states,  state-to-state  vibrational  transfer  rates  are  needed  for  each  energy 
transfer  process.  For  the  model  used  for  these  calculations  it  is  assumed  that  only  single 
quantum  transfers  are  significant,  which  is  a  good  assumption  for  the  relatively  low  levels 
of  vibrational  excitation  expected.  To  maximize  the  use  of  available  vibrational  transfer 
rate  information,  rate  parameterizations  are  used  for  the  intramodal  VV  and  VT 
processes.  The  parameterization  used  for  the  VV  processes  follows  the  method  of 
Kurnosov  et  al[2. 1 .7].  Those  used  for  VT  processes  follows  the  work  of  Dunnwald  et 
al[2.1 .6].  The  parameterizations  are  fits  to  the  available  transfer  rate  data  using  analytical 
expressions  based  on  modified  forms  of  SSH  theory. 

Results  of  the  master  equation  vibrational  transfer  model  give  the  vibrational 
population  as  a  function  of  time.  With  knowledge  of  the  states  emitting  within  the 
collection  bandwidth  of  the  camera,  and  the  Einstein  A  coefficients  for  those  states,  the 
FQY  can  be  obtained  to  the  model  results.  To  aid  in  the  design  of  experiments  the  model 
is  run  for  a  range  of  temperature  and  pressure  conditions  resulting  in  a  two-dimensional 
map  of  FQY. 

CO?  Fluorescence  Pressure  and  Temperature  Dependence 

The  dependence  of  the  CO2  fluorescence  on  pressure  and  temperature  can  be 
determined  by  substituting  the  modeling  results  into  Eqn.  (2.1.1).  To  simplify  this  process 
we  can  rewrite  Eqn.  (2. 1 . 1 )  as  follows 

Sf=C^-cr(P,T,g{y))^{P,T,zl)  (2.1.3) 

where  C  is  a  constant  that  can  be  determined  through  calibration  with  a  known 
concentration  of  CO2.  Equation  (2.1.3)  indicates  that  we  can  elucidate  the  temperature 
and  pressure  dependencies  of  the  fluorescence  signal  by  plotting  the  product  of  the 
absorbing  species  number  density,  the  absorption  cross  section,  and  the  fluorescence 
quantum  yield.  Figure  (2.1.1)  shows  a  plot  of  this  product  as  a  function  of  temperature 
and  pressure  for  pure  CO2.  The  temperature  and  pressure  dependence  of  the  fluorescence 
is  relatively  small  over  a  large  portion  of  the  temperature  and  pressure  range,  indicating 
that  the  excitation  frequency  used  in  the  plot  would  be  a  good  choice  for  imaging  the 
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species  mole  fraction.  We  can  eliminate  the  dependence  on  local  mole  fraction  by  taking 
the  ratio  of  an  offset- wavelength  excited  image  and  a  online-wavelength  excited  image. 


100  200  300  400 

Pressure,  torr 

Figure  2.1.1.  Temperature  and  pressure  dependence  of  fluorescence  from  C02  as  given  by  the  product 
of  the  number  density,  absorption  cross  section,  and  fluorescence  quantum  yield.  Excitation  frequency  of 
4992.516  cm'1. 

Radiative  Trapping 

Infrared  fluorescence  is  always  resonant  with  an  absorption  band  in  the  species  from 
which  the  fluorescence  originated;  therefore,  reabsorption  of  the  fluorescence,  in  the 
mixture  being  excited,  is  always  occurring  to  some  extent.  The  reabsorption  of 
fluorescence  before  it  reaches  the  camera  is  known  as  radiative  trapping.  Trapping  of  the 
fluorescence  has  two  main  influences  on  the  fluorescence  signal,  the  first  effect  being  a 
reduction  in  the  signal  intensity  reaching  the  camera,  and  the  second  effect  being  the 
introduction  of  error  in  the  measurement  due  to  nonuniform  trapping  of  the  fluorescence 
across  the  field  of  view  of  the  camera.  Nonuniform  trapping  results  from  inhomogenieties 
in  temperature,  pressure,  or  species  mole  fraction  between  the  plane  of  fluorescence 
emission  and  the  plane  of  the  imaging  array. 

For  CO2  radiative  trapping  often  has  significant  influence  especially  when  using  4.3 
pm  collection  as  is  used  in  these  experiments.  This  fluorescence  corresponds  to  a  Av  =  1 
change  in  the  V3  vibrational  mode.  This  mode  corresponds  to  the  asymmetric  stretching 
vibration  of  CO2  which  is  one  of  the  strongest  transitions  in  the  infrared  spectral  region. 
Even  for  the  short  distances  used  in  this  experiment,  the  influence  of  radiative  trapping  is 
significant.  The  trapping  of  the  fluorescence  will  be  nonuniform  in  this  instance,  as  the 
edges  of  the  jet  have  lower  CO2  concentration  than  the  core  of  the  jet;  therefore, 
fluorescence  from  the  edges  is  attenuated  less  than  fluorescence  from  the  core  and  the 
trapping  will  serve  to  augment  the  signal  at  the  jet  edges  relative  to  the  core.  For  imaging 
of  CO2  mole  fraction  this  is  problematic;  however,  for  imaging  of  pressure  or  temperature 
using  the  ratio  of  two  images  this  effect  will  cancel,  and  we  can  recover  the  pressure  or 
temperature  without  the  influence  of  nonuniform  trapping. 

Experimental  Setup 

A  small  vacuum  tight  cell  with  a  converging  nozzle  mounted  at  the  top  is  constructed 
to  study  an  underexpanded  CO2  jet  using  IRPLIF.  The  nozzle  for  the  jet  is  a  60° 
converging  nozzle  with  an  exit  diameter  of  1 .5  mm.  A  pulsed  solenoid  valve  is  used  to 
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turn  the  jet  on  and  off.  Images  of  the  CO2  jet  injected  into  N2  are  acquired  shortly  after 
the  pulsed  solenoid  valve  opens.  A  plenum  upstream  of  the  pulsed  valve  is  used  to 
dampen  out  pressure  fluctuation  resulting  from  the  opening  of  the  valve.  Capacitance 
manometers  (Setra  280E)  are  used  to  monitor  the  pressure  in  the  plenum  and  in  the  cell. 
The  cell  is  connected  to  a  gas-handling  manifold  allowing  the  cell  to  be  filled  with 
nitrogen  prior  to  the  CO2  injection  and  the  cell  to  be  evacuated  following  the  injection. 
The  plenum  pressure  for  the  experiments  is  1340  torr.  The  cell  pressure  is  set  to  170  torr. 
This  gives  a  pressure  ratio  (poo/po)  of  0.126  corresponding  to  a  fully  expanded  Mach 
number  of  2.0. 

Laser  excitation  of  CO2  is  provided  by  a  OPO/OPA  system  (Continuum  Mirage  3000) 
that  has  been  modified  with  a  second  set  of  KTA  crystals  in  the  final  amplifier  stage  to 
increase  conversion  efficiency  [2.1.1].  The  laser  outputs  pulses  of  approximately  20  mJ 
with  a  4-6  ns  FWHM  duration  at  the  target  frequency  of  4992.516  cm'1.  As  discussed 
earlier,  the  spectral  width  of  the  laser  is  approximately  500  MHz  FWHM.  The  laser  beam 
is  formed  into  a  sheet  using  a  combination  of  a  500  mm  focal  length  spherical  CaF2  lens 
and  a  250  mm  focal  length  cylindrical  CaF2  lens.  Placing  the  cylindrical  lens  1000  mm 
after  the  focal  point  of  the  spherical  lens  produces  a  slightly  diverging  sheet  height  of 
8.3  mm  FWHM  and  a  400  pm  FWHM  waist. 

Images  are  acquired  with  a  LN2  cooled  InSb  camera  (Santa  Barbara  SBF134).  The 
camera  is  outfitted  with  a  wide  bandpass  (3.0-5.0  pm)  interference  filter  with  a 
transmission  of  approximately  90%  over  most  of  the  transmission  range.  The  filter 
provided  high  rejection  of  the  scattered  laser  light,  as  well  as  rejection  of  some  thermal 
background  signal.  A  25  mm  diameter,  25  mm  focal  length  CaF2  singlet  lens  is  used  to 
focus  the  fluorescence  onto  the  camera  array,  resulting  in  a  magnification  of  0.53.  The 
pixel-by-pixel  gain  and  offset  of  the  camera  with  the  filter  and  lens  are  calibrated  by 
taking  images  of  grey  body  emitters  at  several  temperatures.  The  camera  could  be  gated 
with  integration  times  as  small  as  130  ns,  and  300  ns  integration  times  are  used  for  this 
experiment. 

Imaging  Results 

Images  are  taken  using  both  the  online-excitation  and  offset-excitation  discussed 
earlier.  Figure  2.1.2  shows  a  single-shot  image  of  a  pure  CO2  jet  issuing  into  N2  using 
online-excitation.  The  image  has  been  corrected  for  background  signal  and  the  laser  sheet 
profile  and  has  been  binned  into  two-pixel  by  two-pixel  cells.  Additionally,  a  three-pixel 
x  three-pixel  Gaussian  filter  is  applied  to  smooth  the  data.  As  expected  from  the 
modeling  for  online  excitation,  we  do  not  see  significant  fluctuation  in  signal 
corresponding  to  the  pressure  and  temperature  variations  in  the  flow  field.  Instead,  in  the 
core  of  the  jet,  we  see  a  relatively  uniform  signal  over  most  of  the  image.  Since  online 
excitation  is  relatively  insensitive  to  temperature  and  pressure,  the  signal  reflects 
primarily  changes  in  the  CO2  mole  fraction.  One  notable  exception  is  at  the  very  top  of 
the  image,  closest  to  the  jet  exit  where  the  Riemann  wave  is  located,  where  the  significant 
variation  in  pressure  and  temperature  cause  a  strong  gradient  to  show  up  in  the  image. 
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Figure  2.1.2.  Online  excitation  image  of  a  supersonic  C02  jet  into  N2  with  a  fully  expanded  Mach 
number  of  2.0  issuing  from  a  1 .5  mm  diameter  converging  nozzle. 
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Figure  2.1  J.  Offset  excitation  image  of  a  supersonic  C02  jet  into  N2  with  a  fully  expanded  Mach 

number  of  2.0  issuing  from  a  1.5  mm  diameter  converging  nozzle. 
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An  example  single-shot  image  result  for  offset  excitation  is  shown  in  Fig.  2.1.3.  In  this 
image,  the  pressure  variations  are  readily  apparent,  along  with  the  cellular  nature  of  the 
jet.  The  signal  reflects  both  variations  in  pressure,  as  well  as  C02  mole,  with  some 
dependence  upon  temperature  and  radiative  trapping.  To  cancel  the  signal  dependence 
upon  C02  mole  fraction  and  eliminate  the  trapping  influence,  we  ratio  the  offset- 
excitation  image  by  the  online-excitation  image.  We  can  convert  the  data  to  pressure  if 
we  neglect  the  sensitivity  to  temperature  and  use  an  average  value  of  temperature  instead. 
Assuming  an  average  temperature  of  225  K,  we  obtain  the  pressure  variation  shown  in 
Fig.  2.1.4. 


x/d 

Figure  2.1.4.  Estimated  center  line  variation  in  pressure  using  the  ratio  of  offset  and  online  excitation 
images.  Plot  is  generated  using  an  average  of  ten  images  for  both  the  offset  and  online  images. 


Current  Status 

A  model  for  IRPLIF  of  C02  is  developed  and  applied  to  optimize  an  imaging 
experiment  for  an  underexpanded  jet  with  a  fully  expanded  Mach  number  of  2.0.  Single¬ 
shot  images  of  the  jet  are  obtained  using  a  300  ns  integration  time.  A  technique  is 
developed  to  measure  pressure  from  the  ratio  of  images  of  the  jet  excited  at  line  center  of 
the  selected  rotational  transition  with  the  laser  detuned  by  0.04  cm'1.  The  ratio  of  these 
two  images  gives  the  local  pressure  with  only  small  errors  due  to  the  temperature 
dependence  of  the  fluorescence.  To  our  knowledge,  these  images  represent  the  first 
IRPLIF  images  of  C02  in  a  supersonic  flow.  Future  work  will  focus  on  further 
development  of  the  technique.  For  example,  by  optimizing  the  excitation  wavelengths  it 
should  be  possible  to  reduce  the  temperature  dependence  further,  or,  by  using  the  ratio  of 
two  images  excited  at  different  frequencies,  the  local  static  temperature  in  the  supersonic 
flow  could  be  determined.  Recently,  the  infrared  laser  source  has  been  improved, 
providing  twice  the  excitation  energy. 
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2.2.  Quantitative  LI F  and  PLIF  measurements  of  in  high-  pressure  flames 

Significant  accomplishments  are  made  in  the  development  of  diagnostic  strategies 
and  design  rules  for  quantitative  LIF  and  PLIF  of  NO  in  high  pressure  flames.  In 
addition,  the  first  UV-PLIF  measurements  are  made  on  the  combustion  product  CO2,  and 
this  novel  new  diagnostic  is  demonstrated  in  the  burnt  gases  of  a  high-pressure  methane 
flame. 

2.2.1  Strategies  for  quantitative  NO  LIF 

Nitric  oxide  laser-induced-fluorescence  (NO  LIF)  2-D  imaging  in  high-pressure 
flames  (l-60bar)  is  investigated  using  a  new  multi-spectral  detection  strategy.  This  work 
builds  on  previous  AFOSR-sponsored  research  that  identified  interference  LIF  from  O2 
and  CO2  in  high-pressure  flames  and  optimized  the  choice  of  excitation  strategies  as  a 
function  of  application  conditions.  The  current  work  develops  design  rules  to  optimize 
the  LIF  detection  wavelengths  for  quantitative  2-D  NO-LIF  measurements  over  a  wide 
range  of  pressures  and  temperatures.  Simultaneous  detection  of  LIF  in  multiple 
wavelength  regions  enables  correction  of  the  NO  signal  for  interference  from  O2  and  CO2 
and  allows  simultaneous  imaging  of  all  three  species.  New  experiments  of  wavelength- 
resolved  1-D  LIF  in  slightly  lean  (^=0.9)  and  slightly  rich  (^=1.1)  methane/air  flames  are 
used  to  evaluate  the  design  rules  and  estimate  the  NO  detection  limits  for  a  wide  range  of 
flame  conditions.  The  work  done  here  provides  a  practical  guideline  for  application  of 
instantaneous  1-D  or  2-D  NO-LIF  imaging  strategies  in  high-pressure  combustion 
systems. 

Background 

The  elevated  pressures  of  combustion  environments  in  most  practical  systems 
introduce  significant  complications  in  the  application  of  NO-LIF,  as  identified  a  decade 
ago  [2.2.1].  Recently,  with  AFOSR  support,  we  began  a  comprehensive  study  of 
quantitative  NO-LIF  in  flames  up  to  60  bar,  leading  to  enhanced  understanding  of  high- 
pressure  NO-LIF.  Different  excitation  strategies  are  investigated  and  optimized 
excitation  wavelengths  for  NO-LIF  recommended  for  various  combustion  environments 
in  a  series  of  papers  [2.2.2]. 

During  the  present  grant,  we  investigated  a  new  multi-spectral  detection  strategy  for 
quantitative  2-D  NO-LIF  imaging.  While  the  spectral  analysis  method  is  used  routinely 
in  bio-chemistry  assays  for  complex  polyatomic  molecules,  to  our  knowledge,  this  is  the 
first  application  of  multi-spectral  detection  in  combustion  LIF  and  certainly  the  first  time 
for  NO  PLIF  diagnostics.  This  work  builds  on  the  results  of  our  previous  efforts  and 
leads  to  a  new  strategy  based  on  the  use  of  multiple  detection  wavelength  regions  for 
simultaneous  determination  of  NO  and  the  interference  from  O2  and  CO2  LIF.  Design 
rules  are  developed,  and  the  detection  wavelength  regions  or  “windows”  are  optimized, 
using  new  wavelength-resolved  1-D  NO-LIF  experiments  in  slightly  fuel  lean  (^=0.9) 
and  slightly  fuel  rich  (^=1.1)  methane/air  flames.  Demonstration  measurements  of  2-D 
LIF  imaging  then  are  performed  to  validate  the  imaging  strategy  and  estimate  the  NO 
detection  limits  for  lean  and  rich  high-pressure  flames.  Quantitative  measurements  of  NO 
in  these  seeded  flames  are  compared  with  detailed  chemical  kinetics  calculations  as  a 
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function  of  pressure  and  stoichiometry,  and  good  agreement  is  observed  within  the  limits 
imposed  on  the  calculation  by  the  uncertainty  in  the  measured  temperature. 

Experimental  Details 

The  apparatus  used  in  this  study  has  been  described  previously  [2.2.2]  and  only  is 
discussed  here  briefly.  Laminar,  premixed  methane/air  flat-flames  at  pressures  from 
1  to  60  bar  are  stabilized  on  a  porous,  sintered  stainless  matrix  of  8  mm  diameter.  This 
burner  is  mounted  in  a  pressurized  vessel  with  an  inner  diameter  of  60  mm  with 
stabilization  to  ±0.1  bar.  Investigations  are  conducted  for  fr=0.9  and  1.1  methane/air 
equivalence  ratios.  Optical  access  to  the  flame  is  possible  via  four  quartz  windows 
(Heraeus,  Suprasil  2  Grade).  Laser  pulse  energy  ranged  from  0.8  to  1.2  mJ/pulse  (7  ns 
pulse  at  10  Hz  with  linewidth  of  0.2  cm’1)  from  a  Nd:YAG-pumped  (Quanta  Ray 
GCR250)  frequency-doubled  (BBO)  dye  laser  (LAS,  LDL205).  For  spectrally  resolved 
1  -D  line  imaging  measurements,  the  beam  is  focused  weakly  (diameter  =  1  mm)  along  a 
line  3  mm  above  the  burner  matrix  and  crosses  the  flame  horizontally.  The  pulse  energy 
is  monitored  with  a  fast  photodiode  (LaVision).  Fluorescence  signals  are  collected  at 
right  angles  to  the  laser  beam  and  focused  with  a /  =  105  mm,//4.5  achromatic  UV-lens 
(Nikon),  dispersed  spectrally  through  an  imaging  spectrometer  (Lavision  Chromex 
250IS),  and  imaged  onto  an  intensified  CCD  camera  (LaVision  Dynamight).  For  2-D 
imaging  measurements,  the  laser  is  stretched  optically  into  a  vertical  sheet 
(0.5  mm  x  15  mm)  and  imaged  directly  using  an  ICCD  camera  with  different  filter  sets 
for  suppression  of  Rayleigh  Scattering  and  bandpass  filtering.  The  signals  from  20  laser 
pulses  are  accumulated  on  chip  for  all  measurements.  For  multi-window  imaging 
measurements,  wavelength  regions  are  isolated  with  bandpass  filters  (discussed  below) 
and  separately  collected  on  a  single  camera.  For  practical  2-D  imaging  applications 
where  all  spectral  windows  require  simultaneous  detection,  a  single  laser,  multi-camera 
strategy  can  be  used. 

Measurement  Strategy 

There  are  three  significant  complications  that  arise  when  applying  NO-LIF  to  high- 
pressure  combustion  environments:  (I)  pressure  broadening  and  shifting  of  NO 
transitions,  (2)  strong  attenuation  of  UV  light  and  (3)  interference  LIF.  Excitation 
features  of  NO  broaden  with  increased  pressure,  and  for  pressures  greater  than  10  bar,  it 
is  no  longer  possible  to  excite  single  rotational  features  of  the  excitation  spectrum  due  to 
the  spectral  overlap  [2.2.2].  Along  with  an  increase  in  fluorescence  quenching  of  the  NO 
A  state,  pressure  broadening  leads  to  a  decrease  in  the  overall  signal  strength  since  the 
laser  linewidth  (0.2  cm'1)  is  less  than  the  broadened  linewidths.  Second  is  the  strong 
attenuation  effect  of  both  the  excitation  laser  light  and  fluorescence  signal.  UV  light  at 
short  wavelengths  (<250  nm)  is  absorbed  strongly  in  hot  post-combustion  gases,  due 
mainly  to  strong  absorption  by  CO2  and  H2O  [2.2.3]  (CO2  being  the  more  dominant 
contributor  in  our  wavelength  range  and  temperatures).  This  effect  increases  with 
pressure  as  densities  of  the  absorbing  species  increase.  Finally,  there  can  be  significant 
signal  interference.  The  Schumann-Runge  bands  of  O2  are  pervasive  in  the 

same  excitation  and  detection  region  and  can  be  a  major  interference  source  in  fuel  lean 
and  non-premixed  flames  [2.2.1].  Recently  identified  C02  broadband  emission 
(200-450nm)  is  also  a  concern  in  high-pressure  flames  [2.2.4].  Additionally,  polycyclic 
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aromatic  hydrocarbon  (PAH)  LIF  emission  can  be  a  source  of  significant  interference  in 
fuel-rich  flames. 

Isolation  of  NO-LIF  signal  from  interference  sources,  mainly  from  LIF  of  hot  O2  and 
CO2  is  investigated  using  LIF  detection  in  multiple  wavelength  windows,  a  novel  concept 
in  NO-LIF  imaging.  While  an  ultimate  goal  of  our  work  is  the  development  of  a  single¬ 
shot  2-D  imaging  strategy,  all  the  data  presented  here  are  obtained  by  accumulating 
multiple  laser  pulses  on  an  ICCD  camera  (20  shots/image),  made  possible  by  the  steady 
nature  of  our  test  flame.  By  averaging  over  multiple  laser  pulses,  we  are  able  to  estimate 
the  single-shot  LIF  signals  achievable  using  lasers  above  the  modest  power  range  (0.8  to 
1.2  mJ/pulse)  of  our  current  setup.  Figure  2.2.1  shows  a  typical  emission  spectrum  of 
NO-LIF  using  I-D  imaging  with  ^-X(0,0)  band  excitation,  where  non-linear  fitting  and 
computational  simulation  of  NO,  02,  and  C02  LIF  emission  are  used  to  separate  the 
spectrum  into  components  from  NO,  02,  and  C02  LIF  and  the  Rayleigh  scattering  of  the 
laser  .  It  can  be  seen  that  O2  and  CO2  LIF  are  emitted  throughout  the  entire  region  of 
interest  and  overlap  with  the  spectral  features  of  NO. 

Filter  1  Filter  2  Filter  3 


Figure  2.2.1.  Emission  Spectra  of  NO-LIF  with  226.03nm  Excitation  (40bar  flame,  <jf= 0.9,  300  ppm  NO 
seeding).  Spectra  are  obtained  by  taking  I  -D  line-imaging  (top  box).  NO,  O,  and  C02  emission  profiles 
are  enhanced  for  clarity  (bottom  3  boxes).  Filter  I,  filter  2  and  filter  3  indicate  the  detection  windows. 
Profiles  are  taken  in  a  1mm  region  in  the  center  of  the  flame. 

Spectrally  resolved  1-D  imaging  is  obtained  by  aligning  the  laser  beam  parallel  to  the 
spectrometer  slit,  and  the  LIF  signal  is  dispersed  onto  an  ICCD  camera.  This  camera 
image  has  the  position  along  the  laser  beam  on  the  ordinate  and  the  fluorescence 
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wavelength  on  the  abscissa.  1-D  line-imaging  provides  increased  signal-to-noise  ratio 
(SNR)  ratio  due  to  higher  laser  intensity  (beam  vs.  sheet),  direct  visualization  of 
interference  sources,  and  a  more  flexible  selection  of  excitation/detection  strategies,  since 
NO-LIF  can  be  isolated  easily  from  interferences.  In  addition,  only  a  single  laser  and 
single  camera  detection  setup  is  required. 

With  modern-day  high-power,  short  pulse  (<10ns)  lasers,  LIF  signal  dependence  on 
laser  energy  can  show  non-linear  effects  due  to  saturation  of  the  excitation  transition.  In 
our  1-D  line  imaging  experiments  we  observe  only  a  slight  deviation  from  linear 
excitation  when  the  laser  spectral  intensity  (7V)  reach  ~30  W/mV  (I  mJ,  7  ns  pulse  @ 
226.03  nm  with  0.2  cm'1  bandwidth  and  a  focused  diameter  of  1  mm)  in  an  atmospheric 
pressure  flame.  As  pressure  is  increased  to  the  pressures  of  current  interest,  the  rapid 
increase  in  the  rotational  energy  transfer  refills  the  depleted  ground  state  while  removing 
excess  population  in  the  upper  excited  state,  thereby  greatly  extending  the  linear  regime 
of  NO-LIF.  This  population  recycling  effect  enables  the  use  of  higher  laser  power  for 
linear  single-shot  imaging  in  high-pressure  environments  (e.g.,  up  to  20  mJ/pulse  @ 
20  bar  for  our  setup)  resulting  in  increased  SNR  ratios.  For  2-D  imaging  measurements, 
significantly  stronger  laser  pulses  can  be  used  since  the  laser  beam  is  dispersed  into  a 
sheet  and  the  spectral  intensity  is  decreased. 

Experimental  Results 

Figure  2.2.1  shows  the  detection  windows  (filters)  used  for  the  2-D  imaging 
measurements.  Filter  I  (254  nm,  FWHM  17.5  nm;  Asahi  Optics)  is  the  primary  detection 
channel  for  NO,  chosen  to  optimize  the  NO  LIF  signal,  while  filter  2  (280  nm,  FWHM  12 
nm;  Lambda  Research  Optics)  and  filter  3  (340  nm,  FWHM  10  nm;  Asahi  Optics)  are 
used  to  correct  the  signal  for  interference  from  O2  and  CO2,  respectively.  In  addition,  a 
240  nm  long-pass  filter  is  used  for  additional  suppression  of  the  Rayleigh  scattering.  The 
selection  rules  for  filter  1,  the  primary  NO-LIF  detection  window,  maximize  the 
selectivity  and  intensity  while  minimizing  the  contribution  from  Rayleigh  scattering  of 
the  laser.  The  bandpass  width  also  should  be  wide  enough  to  avoid  bias  in  the  collection 
efficiency,  depending  on  the  rotational  levels  excited.  For  our  test  flame,  optimization  is 
carried  out  by  taking  the  spectrally  resolved  line-imaged  data  (Fig.  2.2.1)  and  scanning  a 
simulated  Gaussian  filter  function  (with  varying  spectral  widths)  throughout  the 
fluorescence  spectrum.  For  detection  windows  (filter  2  and  filter  3)  that  are  used  to 
isolate  interference  LIF,  the  selection  criteria  require  regions  with  maximum  fractional 
contribution  from  the  specific  interference  species  for  increased  accuracy  in  the 
correction. 

Selection  of  the  optimum  number  of  detection  wavelength  windows  requires 
knowledge  regarding  the  global  chemistry  of  the  combustion  process  involved.  For  many 
experimental  applications,  single-window  detection  (filter  I)  is  sufficient,  as  the  flames 
have  minimal  O2  and  CO2  interference  if  the  proper  excitation  and  fluorescence 
wavelengths  are  selected.  Non-sooting  flames  with  equivalence  ratios  near  one  or  slightly 
richer  generally  exhibit  negligible  levels  of  signals  from  02,  and  contribution  from  C02 
LIF  also  can  be  neglected  for  pressures  less  than  20  bar.  For  high-pressure  flames 
(—20  bar)  with  lean  chemistry  (^<0.9),  the  interference  from  CO2  and  O2  requires  more 
than  one  detection  window  for  sensitive  and  accurate  detection  of  NO.  In  the  case  of  rich 
flames  (<£>1.5)  or  flames  with  incomplete  combustion  of  the  fuel,  interference  from  PAH 
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and  other  hydrocarbons  additionally  can  complicate  NO-LIF  measurements  and  will 
require  further  work  to  establish  an  appropriate  detection  strategy. 

Accurate  2-D  NO-LIF  images  require  corrections  for  various  effects  including  laser 
and  fluorescence  signal  attenuation,  which  vary  with  gas  temperature.  For  steady 
laboratory  flames,  time-consuming  strategies  such  as  multi-line  NO-LIF  fitting  [2.2.5] 
can  be  used  to  obtain  accurate  2-D  temperature  fields.  For  practical  applications,  2-line 
thermometry  offers  a  rapid  means  to  obtain  2-D  temperature  images  with  relatively  high 
precision.  The  temperature  fields  used  in  this  study  are  obtained  via  2-D  2-line 
thermometry  of  NO-LIF  (line  selection:  226.03  nm  and  225.847  nm)  and  illustrated  in 
Fig.  2.2.2.  A  detection  window  centered  at  254  nm  (FWHM  17.5  nm;  Asahi  Optics)  is 
used  to  maximize  the  signal  purity  of  the  NO-LIF. 


Figure  2.2.2.  NO  Temperature  fields  from  2-D  2-line  NO-LIF  thermometry  (line  selection  :  226.03  nm 
and  225.847  nm)  for  various  pressures  (1-60  bar)  at  ^=0.9.  Obtained  by  accumulating  20  shots  on-chip 
and  using  2-D  PLIF  imaging  of  NO-LIF. 

Quantitative  2-D  NO-LIF  images  obtained  from  the  multi-window  detection  strategy 
are  shown  in  Fig.  2.2.3.  The  data  at  pressures  of  1  and  10  bar  use  a  two-window  (filter  I 
&  2)  detection,  since  interference  from  CO2  is  negligible  under  these  conditions  [2.2.4], 
The  second  window  provides  a  means  to  correct  the  raw  image  for  the  interference  of  O2- 
LIF  in  slightly  lean  conditions  (^=0.9).  In  slightly  rich  conditions  (^=1.1),  02  is  mostly 
consumed  during  the  reaction,  and  single-window  (Wl)  detection  is  sufficient  (NO-LIF 
signal  selectivity  >98%)  and  no  correction  is  needed.  For  pressures  of  20  bar  and  higher, 
a  three-window  detection  strategy  is  used  for  correction  of  the  raw  images  for  both  O2 
and  C02  LIF.  The  additional  filters  allow  simultaneous  images  of  hot  CO2  and  02.  The 
data  in  Fig.  2.2.4  show  the  image  of  02  and  C02  extracted  from  the  interference  signals 
for  rich  and  lean  40  bar  flames  (right  side)  along  with  NO  images  with  and  without  these 
corrections  (left  side).  It  should  be  noted  that  the  impact  of  the  correction  on  the  NO 
intensity  is  not  significant  due  to  the  fact  that  the  excitation  wavelength  used  already  has 
been  optimized  for  maximum  suppression  of  the  interference.  The  data  in  Fig.  2.2.3  show 
that  for  slightly  lean  flames  (<j>=0.9),  NO  concentrations  increase  with  pressure  as 
predicted  by  the  thermal  and  prompt  NO  formation  mechanisms.  For  slightly  rich 
conditions  (<|>=I.l),  NO  rebum  due  to  excess  hydrocarbons  in  the  flame  causes  the 
concentrations  to  decrease  with  pressure. 

Current  Status 

Spectrally  resolved  l-D  line-imaging  offers  unique  advantages  in  applying  single¬ 
shot  detection  of  NO-LIF  to  practical  systems.  In  the  current  work,  l-D  line-imaging 
measurements  are  carried  out  for  a  range  of  flame  pressures  (1-60  bar)  and 
stoichiometries  to  provide  guidelines  for  filter  selection  in  the  multi-window  detection 
strategies  used  for  2-D  measurements.  These  l-D  measurements  do  not  require  the  laser 
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beam  to  be  spread  into  a  sheet,  and  the  higher  laser  energy  intensity  results  in  increased 
SNR.  This  work  provides  practical  guidelines  for  application  of  instantaneous  1-D  or  2-D 
NO-LIF  imaging  in  high-pressure  combustion  systems.  Further  studies  to  identify  the 
minimum  detectable  of  NO  concentrations  using  1-D  line-imaging  are  underway. 


Figure  2.23.  Quantitative  NO  concentrations  from  multi-window  2-D  NO-LIF  Imaging  (1-60  bar,  ^=0.9 
and  1.1  with  300  ppm  NO  seeding).  Note  that  the  NO  concentrations  in  rich  flames  are  significantly 
reduced  by  NO  rebum. 
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Figure  2.2.4.  NO,  02  and  C02  concentration  fields  from  a  three-window  2-D  imaging  strategy  applied  in 
a  40  bar  flame  (4>=0.9  and  1,1  with  300  ppm  NO  seeding).  For  NO,  interference  corrected  and  uncorrected 
images  are  both  shown  (left  panel).  LIF  Intensities  are  calibrated  using  Chemkin  Premix  1-D  flame 
calculations. 
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2.2.2  First  PLIF  imaging  of  C02 


UV  planar  laser-induced  fluorescence  (PLIF)  images  of  hot  carbon  dioxide  (CO2)  are 
obtained  in  a  laminar  flame  (CH^air)  at  high  pressure  (20  bar)  with  excitation 
wavelengths  at  239.34  nm  and  242.14  nm.  Excitation  wavelengths  are  chosen  to 
minimize  the  contribution  of  nitric  oxide  and  molecular  oxygen  LIF  signals.  Spectrally 
resolved  single-point  measurements  are  used  for  correction  of  remaining  oxygen  LIF 
interference.  The  continuum  LIF  signal  from  electronically  excited  CO2  is  detected  in  a 
broad  (280-400  nm)  emission  region.  The  UV  PLIF  of  hot  C02  has  the  potential  for 
application  to  a  wide  variety  of  diagnostic  needs  in  high-pressure  flames,  combustors, 
and  engines. 

Background 

At  room  temperature  CO2  is  transparent  for  wavelengths  above  200  nm,  but  recently, 
with  AFOSR  support,  we  discovered  that,  at  flame  temperatures,  the  absorption  of  hot, 
vibrationally  excited  C02  extends  beyond  300  nm  [2.2.3,  2.2.6].  The  importance  of  UV 
absorption  by  hot  C02  for  quantitative  UV  laser  diagnostics  in  combustion  systems  was 
recognized  only  recently  [2.2.3].  Subsequently,  it  was  discovered  that  some  of  the  CO2 
molecules  excited  by  this  absorption  emit  a  broadband  LIF  with  a  faint,  superimposed 
structure  between  200  and  400  nm  [2.2.4].  This  broadband  LIF  emission  was  found  to  be 
present  in  spectrally-resolved  measurements  in  CfVair  and  CHVAr/Ch  flames  at  elevated 
pressures,  where  the  intensity  of  the  LIF  signal  varied  with  CO2  absorbance  and  the  lack 
of  signal  in  a  20  bar  H2/CVAr  flame  eliminated  the  possibility  of  LIF  from  H20.  It  was 
shown  that  the  C02  signal  is  linear  with  number  density  and  laser  fluence  within  the 
investigated  ranges  [2.2.4]. 

During  the  present  AFOSR-grant  period,  we  acquired  the  first  UV  LIF  images  of  C02 
in  a  laminar,  steady  high-pressure  flame  (20  bar).  Although  these  first  images  are  not 
quantitative  concentration  maps  of  CCh,  they  illustrate  the  potential  to  image  the 
distribution  of  hot  exhaust  gases  in  practical  high-pressure  combustors. 

CO?  UV  Planar  Laser  Induced  Imaging 

To  image  C02  in  high-pressure  flames,  a  laminar,  premixed  methane/air  flame  is 
stabilized  on  a  porous  sintered  stainless  steel  matrix  of  8  mm  diameter.  This  facility  was 
developed  with  AFOSR  support  and  now  is  used  for  our  measurements  of  NO-LIF  at 
high-pressure  described  in  section  2.2.1. 

Detection  of  C02  in  high-pressure  flames  following  absorption  of  a  photon  in  the 
200-300  nm  range  is  complicated  by  the  simultaneous  excitation  of  NO  in  the  A2L+-X2n 
y  bands  and  02  in  the  B3ST-X5I+  Schumann-Runge  system.  Although  pressure 
broadening  increases  the  overlap  in  these  spectra,  an  excitation  wavelength  can  be 
selected  that  does  not  overlap  with  any  NO  transitions;  however,  due  to  the  pervasive 
nature  of  multiple  vibrational  bands  of  02  in  this  region,  it  is  impossible  to  avoid  02-LIF 
interference  completely.  We  choose  a  strategy  to  use  excitation  wavelengths  which 
avoids  overlap  with  NO  transitions  in  the  region,  while,  at  the  same  time,  minimizing  the 
overlap  with  the  multiple  vibrational  bands  of  02.  Two  candidate  excitation  wavelengths 
are  selected,  239.34  nm  and  242.14  nm,  both  of  which  are  located  between  the  A-X '(0,1) 
and  A— X  (0,2)  bands  of  NO,  as  shown  in  Fig.  2.2.5.  These  choices  minimize 
contributions  of  02,  as  predicted  by  numerical  simulations  and  confirmed  through 
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experimental  excitation  scans  with  wavelength-resolved  detection  in  a  lean  high-pressure 
flame  ( <f>  =  0.8,  p  =  10  bar),  where  excess  O2  yields  identifiable  L1F  signal  in  the  burnt 
gas.  In  comparison,  excitation  at  242.14  nm  is  found  to  provide  the  least  interference 
from  O2  with  sufficient  NO  minimization,  while  239.34  nm  excitation  is  more  selective  in 
completely  avoiding  any  NO  transition. 

The  detection  window  of  the  CO2-LIF  signal  extends  from  280  to  400  nm.  The 
separation  of  this  spectral  region  from  the  laser  wavelength  ensures  suppression  of 
Rayleigh  scattering,  as  well  as  Raman  scattering  from  other  combustion  products.  The 
400  nm  blue  limit  also  suppresses  background  from  flame  emission,  scattered  pump-laser, 
and/or  room  light.  Fig.  2.2.6  shows  a  spectrally-resolved  CO2  UV-L1F  emission  spectrum 
in  the  center  region  and  outer  edge  of  the  flame.  The  spectral  detection  window  is  shown 
as  a  box  in  Fig.  2.2.6.  The  CO2-LIF  signal  corresponds  to  the  broadband  signal  evident 
throughout  the  entire  detection  region,  although  additional  emission  from  O2LIF  appears 
as  we  move  towards  the  outer  rim  of  the  flame,  where  additional  air  is  entrained  and  O2 
levels  are  increased. 


Excitation  Wavelength  [nm] 

Figure  2.2.5.  Simulated  LIF  excitation  spectra  of  02  (upper  frame)  and  NO  (lower  frame)  at  20  bar  and 
2000  K.  A  magnified  02  spectrum  near  the  excitation  wavelengths  also  is  shown.  239.34  nm  and 
242.14  nm  excitation  lines  are  shown  with  arrows. 

The  resulting  images  are  corrected  for  laser  sheet  inhomogeneities  and  the  signal 
attenuation  along  the  path  of  the  laser  beam  caused  by  hot  CO2  and  H2O  absorption. 
Correction  of  attenuation  utilizes  recently  quantified  absorption  cross-section  data  [2.2.3] 
of  hot  CO2,  as  well  as  flame  temperatures  obtained  from  multi-line  fitting  of  NO-LIF 
excitation  scans  [2.2.7].  Additional  correction  for  O2-LIF  contribution  is  carried  out  by 
taking  spectrally  resolved  images  and  extracting  spatial  distribution  of  O2  across  the 
flame  and  subtracting  it  from  the  CO2-LIF  image.  For  this,  spectrometer  measurements 
are  taken  at  a  location  3  mm  above  the  burner  matrix  using  a  beam  instead  of  an 
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expanded  laser  sheet.  The  slit  of  the  spectrometer  provides  the  spatial  resolution  across 
the  horizontal  span  of  the  flame.  Spatial  distribution  of  CO2  across  the  flame  can  be 
extracted  from  the  spectrally  resolved  images  by  assuming  that  the  signal  from  the  center 
region  of  the  flame  is  mainly  CO2.  As  we  move  towards  the  edge  of  the  flame,  additional 
O2-LIF  signal  is  observed  and  can  be  isolated  to  yield  the  amount  of  O2-LIF  interference. 
The  O2-LIF  field  is  extrapolated  to  the  remainder  of  the  flame  using  geometrical 
arguments  and  assuming  uniform  composition  of  the  burnt  gas. 


Figure  2.2.6.  Spectral  detection  bandpass  (280-400  nm)  shown  on  spectrally-resolved  C02  UV-LIF 
spectra  with  242.14  nm  excitation. 

Results 

CO2  visualization  using  UV-LIF  in  a  20  bar  flame  along  with  signal  attenuation  along 
the  path  of  the  laser  and  O2  contribution  information  are  presented  in  Fig.  2.2.7.  Signal 
attenuation  across  the  centerline  can  amount  to  over  10%  at  20  bar  across  a  flame  with  an 
exit  diameter  of  8  mm,  as  illustrated  in  image  a  of  Fig.  2.2.7  [2.2.3].  Although  not 
required  for  our  small  flame,  absorption  of  fluorescence  signal  also  may  need  to  be 
considered  as  it  radiates  from  the  plane  of  the  excitation  light  sheet  to  the  collecting  lens. 
The  attenuation  correction  requires  knowledge  of  the  temperature  field,  which  is 
determined  for  this  flame  using  a  multi-line  NO-LIF  temperature  imaging  technique 
[2.2.5].  For  the  work  reported  here,  attenuation  corrections  are  made  for  the  excitation 
laser  beam  using  Beer-Lamberf  s  law  with  absorption  cross-section  data.  A  pixel  by  pixel 
attenuation  correction  across  the  centerline  of  the  flame  is  made,  and  image  a  of  Fig. 
2.2.7  shows  the  magnitude  of  attenuation  calculations. 

Image  b  is  the  interference  O2-LIF  for  excitation  at  239.34  nm,  and  data  for  242.14 
nm  excitation  are  similar.  The  observed  O2-LIF  contribution  in  the  stable  center  region 
of  the  flame  is  around  1%  of  the  entire  LIF  signal  for  our  flame  conditions  (equivalence 
ratio  of  <j>  =  1.1),  which  we  attribute  to  noise  in  the  extraction  process.  The  O2 
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concentration  increases  around  the  edges  of  the  flame,  where  air  from  the  secondary  flow 
is  entrained  into  the  system. 


Figure  2.2.7,  Results  of  imaging  analysis  of  C02  UV-PLIF  a  total  attenuation  due  to  hot  C02  and  H20 
(239.34  nm  excitation);  b  02-LIF  contribution;  c  and  d  C02-L1F  and  02-LIF  combined  image  for  239.34 
nm  and  242.14  nm  excitation,  respectively;  e  and  f  final  C02-LIF  image  for  239.34  nm  and  242.14  nm 
excitation,  respectively  with  02-LIF  correction.  Images  b-f  share  same  scale  and  image  a  is  ten  times 
more  sensitive. 


Images  e  and  / of  Fig.  2.2.7  show  the  final  CO2-LIF  distribution  using  UV  LIF  with 
excitation  at  239.34  nm  and  242.14  nm,  respectively.  CO2-LIF  intensities  remain  fairly 
constant  to  a  first-order  estimate  throughout  the  hot  burnt  gas  region.  Previous  studies 
with  spectrally  resolved  point  measurements  show  that  the  relative  contribution  of  CO2- 
LIF  signals  increases  with  pressure  and  can  become  a  significant  interference  source  in 
the  detection  of  other  minor  species  in  post  combustion  gases. 

Current  Status 

We  present  2D-L1F  images  of  hot  CO2  using  UV  excitation  in  the  burnt  gases  of  a 
premixed,  high-pressure  laminar  flame.  To  our  knowledge,  these  are  the  first  CO2  images 
published  using  UV-excitation  and  detection.  Excitation  at  239.34  nm  and  242.14  nm 
effectively  can  avoid  overlap  with  NO  transitions  in  the  region  and  minimizes  any 
overlap  with  the  multiple  vibrational  bands  of  O2.  A  first-order  correction  for  residual 
interference  from  O2  LIF  has  been  performed  by  using  spectrally  resolved  emission 
measurements  of  the  flame.  Further  understanding  of  spectral  features  exhibited  by  CO2 
will  be  required  to  enable  a  more  detailed  and  quantified  isolation  of  the  O2  contribution. 

The  CO2  excitation  is  strongly  temperature-dependent,  and  this  visualization  strategy 
will  isolate  the  high  temperature  regions  of  the  combustion  effluent,  possibly  proving  a 
measurement  strategy  to  infer  temperature  in  flames  with  uniform  CO2  concentrations. 
The  strategy  reported  here  for  UV  visualization  of  high-temperature  CO2  in  high-pressure 
combustion  exhaust  also  may  prove  very  useful  to  understand  mixing  phenomena  in  a 
variety  of  practical  combustion  devices. 
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Further  work  is  planned  to  image  CO2  at  higher-pressures  (~60  bar)  and  to  quantify 

the  temperature  and  wavelength-dependence  of  the  CO2  fluorescence  quantum  yield. 

Correction  methods  utilized  in  the  current  study  are  limited  to  stable  flames,  and 

additional  investigation  into  single-shot  measurement  techniques  will  also  be  addressed. 
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2.3.  Temperature  diagnostics  based  on  CO 2  absorption  of  UV  light 

A  diagnostic  for  microsecond-time-resolved  temperature  measurements  using 
ultraviolet  laser  absorption  of  vibrationally  hot  C02  is  developed.  The  absorption  ratio  of 
continuous  wave  laser  radiation  near  244  and  266  nm  is  used  to  determine  temperature. 
Precise  temperature-dependent  absorption  cross  sections  are  measured  for  C02  at  the  two 
laser  wavelengths.  Measurements  behind  shock  waves  in  both  non-reacting  and  reacting 
(ignition)  systems  are  made  and  comparisons  to  isentropic  and  constant  volume 
calculations  are  performed  to  validate  the  temperature  measurement  technique. 

Background 

At  room  temperature  C02  is  transparent  at  wavelengths  longer  than  205  nm,  but  at 
high  temperature  (>1000  K)  C02  has  significant  absorption  in  the  190-310  nm  spectral 
range  due  to  vibrationally  excited  molecules.  With  AFOSR  support,  we  have  investigated 
the  feasibility  of  a  thermometry  technique  using  the  strongly  temperature-dependent  C02 
absorption  cross-section  in  this  region  [2.3.1].  We  first  measured  the  ultraviolet 
absorption  cross-section  of  C02  in  the  wavelength  region  216.5-306  nm  [2.3.2].  The 
spectrally  smooth  nature  of  the  absorption  and  its  strong  temperature-dependence  suggest 
the  use  of  two- wave  length  absorption-based  thermometry  for  high-temperature  gases 
where  C02  is  present.  Additionally,  the  absorption  cross-section  shows  no  pressure 
dependence  in  experiments  up  to  8  atm,  and  none  is  expected  at  higher  pressures  due  to 
its  smooth  pre-dissociative  nature,  enabling  high-pressure  temperature  measurements. 

The  C02  absorption  cross-section  a(X,7)  is  independent  of  pressure  in  the  wavelength 
range  of  interest  and  can  be  described  using  a  semi-empirical  form 


In  o(X,  T)  =  a  +  bX 


(2.3.1) 


where  a  =  c\  +  c2T  +  C3/T  and  b  =  d\  +  d2T  +  dilT.  The  CC^  cross-section,  o(X,  7),  is  in 
units  of  10'19  cm2molecule  ',  the  wavelength,  X,  is  in  units  of  100  nm,  and  the 
temperature,  T,  is  given  in  units  of  1000  K.  The  following  parameters  from  Oehlschlaeger 
et  al.[2.3.2]  provide  a  best  fit  to  the  absorption  cross-section  measurements:  ci  =  0.05449, 
c2  =  0.13766,  o  =  23.529,  d\  =  1.991,  d2  =  -0.17125,  and  d3  =  -14.694.  This  fit  can  be 
used  over  the  range  of  the  cross-section  measurements:  1 500-4500  K  and  2 1 6.5-306  nm. 

Taking  the  ratio  of  absorbance  at  two  different  laser  wavelengths  gives 


absorbance^ t,)  =  <t(/1„7>C02Z  =  c7(A,,r)  = 
absorbance (z^)  a(A7,T)nC02L  0(^,7) 


(2.3.2) 


• 

allowing  the  measurement  of  temperature  independent  of  knowledge  of  C02 
concentration  or  path  length.  For  the  measurements  presented  here  we  chose  to  use  cw- 
laser  radiation  at  244  and  266  nm.  These  two  wavelengths  can  be  generated  conveniently 
and  simultaneously  with  available  laser  sources,  avoid  most  major  interference 
absorption,  and  have  an  appropriate  wavelength  separation  for  measurement  of 
temperature  in  the  region  of  interest  in  shock  wave  experiments  and  combustion  studies. 
A  polynomial  expression  for  temperature  as  a  function  of  the  absorbance  ratio  (cross- 
section  ratio)  from  1500-4500  K  is  given  by 

T=  (2.857  x  1 04)/?4  -  (3.043  x  104)#3  +  (1.478  x  104)/?2  +  (9.57  x  102)/?+  1041  [K] 

(2.3.3) 
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where  R  is  the  ratio  of  266  to  244  nm  absorbance  ( R  =  absorbance(266  nm)  /  absorbance 
(244  nm)). 

The  CO2  absorption  cross-section  at  244  and  266  nm  is  plotted  versus  temperature  in 
Fig.  2.3.1.  The  uncertainty  in  temperature  measurements  made  using  the  absorbance  ratio 
method  is  limited  by  the  magnitude  of  the  fractional  absorption,  particularly  at  the  more 
weakly  absorbing  wavelength  266  nm.  The  measured  ratio  has  greater  uncertainty  at  low 
temperature  due  to  smaller  absorbance  (lower  SNR);  therefore,  the  inferred  temperature 
is  more  uncertain  at  lower  temperatures.  An  uncertainty  in  absorbance  due  to  SNR 
limitations  of  ±0.1%  (-ln(///o)  =  ±0.001)  generally  is  achievable  in  shock  wave 
experiments  regardless  of  the  magnitude  of  the  absorbance.  The  uncertainty  in  measured 
temperature  is  mapped  as  a  function  of  temperature  and  the  product  of  partial  CO2 
pressure  and  path  length  (PcoiL)  in  Fig.  2.3.2.  The  uncertainty  decreases  at  higher 
temperature  and  pressure.  The  uncertainties  given  in  Fig.  2.3.2  assume  an  uncertainty  of 
±0.1%  in  absorbance  (-ln(///o)  =  ±0.001),  and  the  contours  represent  the  percentage 
uncertainty  in  the  measured  temperature.  The  favorable  uncertainty  contours  in  Fig.  2.3.2 
demonstrate  the  attractiveness  of  this  diagnostic  for  high-pressure,  high-temperature 
measurements. 


Figure  2  J.l.  Measured  C02  absorption  cross-section  at  244  and  266  nm. 

Demonstration  of  Time-Resolved  Temperature  Measurements 

The  technique  to  measure  time-resolved  temperature  using  the  ratio  of  CO2 
absorption  at  two  UV  wavelengths  is  demonstrated  behind  reflected  shock  waves  in  a 
pressure-driven  (helium),  high-purity  stainless  steel  shock  tube,  as  illustrated  in  Fig. 
2.3.3.  The  incident  and  reflected  shock  conditions  are  calculated  using  the  ideal  normal 
shock  equations.  Uncertainties  in  the  experimental  pressure  and  temperature  given  by  the 
shock  equations  are  <1%  in  the  CCVAr  mixtures  and  <2%  in  the  CHVCVCCVAr 
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mixtures  used  in  the  reacting  experiments.  Vibrational  equilibrium  is  assumed  behind 
both  the  incident  and  reflected  shock  waves. 


Temperature  (K) 

Figure  23.2.  Uncertainty  in  measured  temperature  for  given  temperature  and  product  of  partial  C02 
pressure  and  path  length  ( PCo2 /•)•  The  contours  represent  lines  of  constant  uncertainty  in  the 
temperature  measurement  assuming  an  uncertainty  in  measured  absorbance  of  ±0.1%  (-ln(/  /  /0)  = 
±0.001). 

The  two-color  UV  absorption  arrangement  is  illustrated  in  Fig.  2.3.3.  Laser  radiation 
at  244.061  nm  (~5  GHz  linewidth)  and  266.075  nm  (<5  MHz  linewidth)  is  produced  by 
the  single  pass  of  a  focused  laser  beam  at  488  (488.122)  nm  (Ar+  line)  and 


■  D 


Figure  233  Schematic  of  experimental  setup:  L,  lens;  BBO,  beta  barium  borate  frequency  doubling 
crystal;  PB,  Pellin-Broca  prism;  M,  mirror;  BS,  beam  splitter;  D,  detector;  PT,  piezoelectric  pressure 
transducer;  and  Kistler,  shielded  Kister  piezoelectric  pressure  transducer. 

532  (532.15)  nm  (Nd:YV04)  through  an  angle-tuned  BBO.  The  harmonics  (244  and  266 
nm)  are  separated  from  the  fundamental  beams  (488  and  532  nm)  with  a  set  of  Pellin- 
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Broca  prisms.  The  two  UV  laser  beams  (1.5  mW  each)  each  are  split  into  two 
components:  one  component,  less  than  I  mm  in  diameter,  passing  through  the  shock  tube 
windows  (UV  fused  silica)  at  a  location  2  cm  from  the  endwall,  and  one  detected  prior  to 
absorption  as  a  reference.  The  beams  are  detected  using  amplified  SI 722-02  Hamamatsu 
silicon  photodiodes  (risetime  <  1.0  ps,  4.1  mm  diameter)  and  recorded  on  a  digital 
oscilloscope. 

Experiments  are  performed  with  several  test  gas  mixtures  to  demonstrate  the 
diagnostic  and  to  characterize  the  effect  of  energy  release  behind  shock  waves. 
Experiments  first  are  performed  in  test  mixtures  of  5  and  10%  CO2  dilute  in  argon. 
Finally,  experiments  are  performed  in  several  stoichiometric  mixtures  of  CHVCVCCh/Ar 
to  demonstrate  the  diagnostic  in  experiments  with  varying  levels  of  energy  release  due  to 
chemical  reactions. 

An  example  non-reacting  experiment  is  shown  in  Fig.  2.3.4  for  a  10%  CCVAr  test 
mixture  (all  data  are  filtered  to  1  MHz  time  resolution).  The  temperature  is  inferred  with 
the  266/244  absorbance  ratio  method  described  above  and  is  found  to  be  nearly  constant 
for  2  ms  before  several  abrupt  changes  produced  by  gas  dynamics.  These  changes  are  due 
to  interaction  of  the  reflected  shock  wave  with  the  contact  surface  (helium/test  gas 
interface)  and  interaction  of  the  reflected  shock  wave  with  the  rarefaction  wave 
(expansion  of  helium  driver  gas).  Additionally,  a  comparison  of  the  measured 
temperature  is  made  with  the  temperature  calculated  from  the  pressure  trace  via  the 
assumption  of  an  isentropic  relationship  between  pressure  and  temperature 


T  = 


(2.3.5) 


Figure  23.4.  Example  non-reacting  experiment:  with  initial  reflected  shock  conditions  from  ideal 
shock  relations:  10%  C02/Ar,  2138  K  and  0.121  MPa.  Measured  temperature  (dark  line)  vs  isentropic 
temperature  (light  line),  measured  temperature  determined  from  absorbance  ratio,  and  isentropic 
temperature  determined  from  pressure  measurement. 
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A  comparison  of  the  temperature  measured  using  the  absorbance-ratio  method  with 
that  resulting  from  the  isentropic  assumption  shows  excellent  agreement  (Fig.  2.3.4),  with 
a  difference  in  the  two  of  no  greater  that  2%  for  this  experiment.  This  good  agreement 
serves  both  to  validate  the  C02  temperature  sensor  and  to  confirm  that  for  shock  wave 
experiments  with  no  energy  release  the  isentropic  assumption  holds. 

Forty-two  experiments  are  conducted  over  a  wide  range  of  temperature  and  CO2 
concentration.  The  scatter  in  the  difference  between  the  measured  temperature  and  the 
ideal  temperature  has  a  I -a  standard  deviation  of  0.6%  and  is  no  greater  than  2.3%.  The 
CO2  absorbance  at  266  and  244  nm  ranges  from  0.5  to  1 3%  and  2  to  32%,  respectively, 
illustrating  the  large  dynamic  range  of  conditions  investigated. 

To  demonstrate  the  C02  absorption  technique  for  thermometry  in  combustion  systems 
and  to  characterize  the  effects  of  energy  release  on  temperature  in  shock  tube  studies,  we 
perform  a  set  of  ignition  experiments  in  CFLt/OyCOyAr  mixtures;  the  time-resolved 
temperature  is  shown  in  Fig.  2.3.5.  Three  mixtures  are  used:  (1)1%  CH4/  2%  02/  10% 
C(V  dilute  in  argon,  (2)  2%  CH4/  4%  02/  10%  CO 2/  dilute  in  argon,  and  (3)  4%  CHJ  8% 
02/  1 0%  CO 2/  dilute  in  argon. 


Figure  2.3.5.  Measured  temperature  (solid  thick  lines)  versus  constant  volume  calculation  (solid  thin 
lines).  Top  graph  (initial  reflected  shock  conditions):  1%  CH4/  2%  02/ 10%  C02/  Ar,  2085  K,  and  0.128 
MPa.  Middle  graph:  2%  CWJ  4%  O -J  1 0%  C02/  Ar,  1 890  K,  0. 1 34  MPa.  Bottom  graph:  4%  CH4/  8% 
O-J  10%  C02/  Ar,  1851  K,  0.126  MPa.  Post-shock  vibrational  equilibrium  is  assumed.  Error  bars 
represent  a  I -a  confidence  interval. 
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These  C02  combustion  measurements  required  experiments  to  investigate  potential 
interference  absorption  from  hot  CH4  and  02.  Absorption  by  CH4  at  these  wavelengths  is 
undetectable,  but  absorption  by  hot  02  is  detectable  and  requires  a  small  correction  to  the 
measurements  at  244.061  nm.  A  chemical  model  and  the  estimated  gas  temperature  are 
sufficient  to  determine  the  small  correction  of  the  absorption  by  02  at  244.061  nm.  This 
interference  can  be  avoided  with  a  different  UV  cw  laser  source  that  can  be  tuned  in 
wavelength  to  an  02  absorption  minimum.  Simulations  of  the  02  spectrum  at  2000  K  and 
0. 1  MPa  in  this  wavelength  range  find  that  the  interference  could  be  reduced  by  a  factor 
of  seven  by  tuning  between  two  absorption  lines  to  a  value  near  243.99  nm. 

The  resulting  temperature  time  histories  are  given  in  Fig.  2.3.5  for  the  three 
combusting  gas  mixtures  along  with  comparison  to  model  results  calculated  using  a 
constant  volume  constraint  and  the  GRI  3.0  methane  combustion  mechanism.  A 
comparison  between  constant  volume  model  calculations  and  measured  temperature 
shows  good  agreement  for  the  ignition  delay  and  the  overall  trends  in  temperature; 
however,  the  post-ignition  temperature  plateau  measured  with  laser  absorption  is  lower, 
in  all  three  cases  than  that  predicted  with  the  constant-volume  model.  This  shortcoming 
of  the  constant  volume  model  is  expected  because  as  energy  is  released  due  to  ignition 
and  the  test  gas  temperature  begins  to  rise  at  the  test  section  location  (2  cm  from  endwall) 
the  hot  combustion  gases  will  expand  down  the  shock  tube,  resulting  in  a  somewhat 
lower  temperature  plateau.  (The  measured  lower  temperature  in  the  post-ignition  plateau 
is  accounted  for  when  determining  the  Q2  interference  absorption  cross-section  at  244 
nm.)  For  the  examples  given  in  Fig.  2.3.5  the  measured  post-ignition  temperature  rise  is 
60-75%  of  that  predicted  by  the  constant  volume  model. 

Current  Status 

We  currently  are  working  to  modify  the  constant  volume  model  and  have  found  good 
success  using  the  CHEMKIN  package  coupled  with  the  measured  pressure  history. 
Initial  experiments  have  yielded  excellent  agreement  with  measured  temperature  time- 
histories.  We  believe  that  this  development  of  a  CHEMKIN  solver  for  reactive  flows 
behind  reflected  shock  waves  is  a  significant  development  that  will  increase  substantially 
the  utility  of  shock  tubes  for  studying  combustion  chemistry. 
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2.4  Advanced  wavelength-multiplexed  diode  laser  absorption  diagnostics 

Results  from  three  different  diode  laser  absorption  sensor  strategies  are  highlighted. 
First,  we  investigate  the  use  of  multiple  absorption  transitions  with  different  temperature 
dependence  to  unravel  non-uniform  temperature  dependence.  The  development  of  an 
absorption  sensing  strategy  for  non-uniform  temperature  distributions  along  the  line-of- 
sight  (LOS)  may  increase  the  potential  applications  of  this  already  successful  sensing 
strategy.  Second,  we  show  the  improved  SNR  for  absorption  sensing  in  the  model 
scramjet  facility  at  Wright  Patterson  Air  Force  Base  using  a  scanned  wavelength- 
modulation  spectroscopy  (WMS)  strategy.  The  results  show  the  importance  of  careful 
optical  design  engineering  and  illustrate  the  potential  of  this  sensing  strategy  for  routine 
use  in  hypersonic  ground  test.  Third,  we  discuss  recent  fuel-sensing  results  using  a  novel 
mid-IR  laser  source  based  on  difference  frequency  generation  (DFG).  This  new  laser  and 
our  differential  absorption  strategy  have  the  potential  for  sensitive  detection  of  a  wide 
variety  of  hydrocarbon  fuels. 

2.4.1.  Non-uniform  temperature  distributions 

Two  measurement  strategies  to  determine  non-uniform  temperature  distributions  in 
combustion  gases  using  line-of-sight  (LOS)  absorption  spectroscopy  have  been  examined 
in  detail.  These  strategies  rely  on  measurements  of  multiple  absorption  transitions  of  a 
single  species,  each  with  unique  temperature  dependence.  The  first  strategy,  called  profile 
fitting,  mathematically  fits  the  observed  absorption  measurements  constrained  with  a 
postulated  temperature  distribution.  The  second  strategy,  called  temperature  binning, 
determines  the  temperature  probability  density  function  (PDF)  along  the  LOS  using 
prescribed  temperature  bins.  The  wavelength-multiplexed  (WM)  sensor  concepts  and  the 
mathematical  representations  are  first  explored  by  simulation  calculations.  Then  these 
concepts  are  demonstrated  in  a  laboratory  arrangement  with  two  distinct  temperature 
regions.  The  measured  “2 -Zone”  T/Xh2o  distribution  is  composed  of  a  hot  flame  zone  (T 
~  1500  K,  Xh20  =  10%)  and  a  cold  room-air  zone  (T  ~  300  K,  Xh20  ~  2%)  each  with  the 
same  pathlength.  The  experimental  results  demonstrate  that  a  non-uniform  temperature 
distribution  can  be  characterized  with  either  strategy.  The  measurement  accuracy  will 
increase  with  the  number  of  transitions  and  also  with  the  use  of  optimally  chosen 
transitions.  The  experimental  results  also  confirm  that  use  of  known  physical  constraints 
to  reduce  the  number  of  degrees  of  freedom  improves  the  interpretation  of  the 
measurement  results  and  thus  the  sensor  performance. 

Background 

Laser  absorption  spectroscopy  has  become  a  successful  sensor  strategy  for  fast,  non- 
intrusive,  sensitive,  and  reliable  in  situ  measurements  of  multiple  flow  field  parameters 
such  as  temperature,  pressure,  velocity,  and  species  concentration  in  a  variety  of  practical 
combustion  and  propulsion  flows  [2.4.1];  however,  this  line-of-sight  (LOS)  sensing 
technique  traditionally  has  been  limited  to  flow  fields  with  nearly  uniform  properties.  In 
practical  flow  fields,  significant  temperature  and  species  concentration  gradients  may 
exist  along  the  measurement  path  for  optical  diagnostics  due  to  flow  mixing,  chemical 
reaction,  phase  change,  heat  transfer  with  the  side  walls,  and  other  effects.  Tomographic 
reconstruction  of  laser  absorption  along  multiple  lines-of-sight  has  been  demonstrated  in 


30 


laboratory  experiments  to  resolve  these  non-uniformities,  but  practical  systems  seldom 
have  sufficient  optical  access  and  the  sensor  redundancy  necessary  for  quantitative 
application  of  tomographic  techniques. 

Extension  of  LOS  absorption  measurements  to  non-uniform  flow  fields  has  been 
explored  previously,  including  work  to  correct  for  boundary  layer  effects,  to  reduce 
sensitivities  to  flow  non-uniformities,  and  to  correlate  pattern  factors  with  the  path- 
averaged  temperatures  inferred  from  different  line  pairs.  Current  efforts  build  on  our 
earlier  AFOSR-sponsored  work  that  first  demonstrated  the  potential  to  extract  the 
temperature  distribution  information  from  simultaneous  measurements  of  multiple 
absorption  transitions  [2.4.2].  Our  current  work  uses  simulation  calculations  to 
investigate  systematically  the  two  strategies:  (1)  profile  fitting  and  (2)  temperature 
binning.  Then  a  demonstration  experiment  is  performed  using  a  wavelength-multiplexed 
sensor  with  an  optical  path  length  that  has  two  distinct  temperature  regions  (~1500K  and 
-300K).  This  experiment  exploits  a  Zolo  fiber-coupled  multiplexer  to  combine  five 
laser  beams  onto  a  single  beam  path  and  a  Zolo  fiber-coupled  de-multiplexer  to  disperse 
the  collected  transmitted  light  onto  five  separate  detectors.  These  five  lasers  interrogate 
seven  absorption  transitions,  which  are  used  to  investigate  a  “2-Zone”  temperature/mole- 
fraction  distribution.  These  results  demonstrate  the  sensor  concept  and  provide  a 
preliminary  investigation  of  sensor  performance. 

Theoretical  Principles 

Building  on  the  ideas  first  demonstrated  in  our  previous  AFOSR-sponsored  work,  we 
first  assume  that  m  (>2)  transitions  have  been  selected  and  their  integrated  Beer’s  Law 
absorbances,  A,,  have  been  determined  from  the  WM-LOS  absorption  measurements.  We 
use  these  simulated  data  to  compare  the  results  from  the  profile-fitting  and  temperature- 
binning  strategies. 

Principles  of  Profile  Fitting 

The  profile-fitting  strategy  first  requires  a  postulated  distribution  of  temperature  along 
the  measurement  path  to  constrain  the  temperature  profile  fitting.  For  example,  in  a 
confined  combustion  flow,  the  gas  temperature  often  has  a  cold  boundary  layer,  and  the 
simplest  representation  would  be  a  “2-T”  profile  with  a  core  flow  at  an  average 
temperature  of  Tc  and  a  boundary  layer  with  an  average  temperature  of  7*  and  a  thickness 
of  Lb,  as  shown  in  panel  a  of  Fig.  2.4.1.  A  more  complex,  but  common,  representation 
would  be  a  parabolic  profile  constrained  by  the  center  temperature  of  Tc  and  the  wall 
temperature  of  Tw,  as  shown  by  panel  b.  Another  somewhat  more  sophisticated  model 
would  be  a  uniform  core  flow  at  a  temperature  of  Tc  and  a  boundary  layer  with  parabolic 
temperature  distribution  constrained  by  the  wall  temperature  of  Tw  and  the  boundary  layer 
thickness  of  Lb,  as  shown  in  panel  c.  These  postulated  temperature  distribution  profiles 
can  be  represented  by  a  general  functional  form  as  follows 

nx)  =  ATchar,L,chor,x),  (2.4.1) 

where  TChar  is  the  characteristic  temperature  such  as  Tc,  Tw  and  7*,  and  Ljchar  is  the 
characteristic  length  such  as  Lb.  Similarly,  the  shape  for  the  absorber  mole  fraction 
distribution  also  is  postulated  using  physical  constraints 

Xabs(X)  =  char>L  Zebari X)  ■>  (2.4.2) 
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where  XChar  is  the  characteristic  mole  fraction  and  Lxchar  the  corresponding  length. 


x  [cm] 


Figure  2.4.1.  Postulated  temperature  distribution  profiles  for  confined  combustion  gases  with  cold 
walls. 

Based  on  the  presumed  temperature  and  mole  fraction  profiles,  the  integrated 
absorbance  of  any  transition  can  be  calculated.  For  the  selected  m  transitions,  a  nonlinear 
equation  set  can  be  established 


A\  ^  8(A ' char'  ^  Xchar’  ^Tchar  ’  x))dx 

^2  =  8(XchariLXchar,X)- S2(f(Tchar,Lrchar,X))dx  ^ 

Am  =  p£g(Xchor, 

L  Xchar  ’  'L Tchar »  xj)dx 


Once  the  number  of  equations,  i.e.  the  number  of  measured  absorption  transitions  m, 
is  larger  than  the  number  of  unknowns,  Eqn.  set  (2.4.3)  can  be  solved  by  nonlinear  least- 
square  fitting 

min  £(  P[g(Xchar,LXchar,x)SXf(Tchar 

’  char  ’  xj)dx-AS  (2.4.4) 

Xchar  ;  =  j  \  / 

to  obtain  the  set  of  TChar,  L  Tchar,  Xchar  and  Lxchar  that  best  describe  the  postulated 
temperature  and  mole  fraction  profiles. 
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Expression  2.4.4  presents  the  most  general  mathematical  model  for  the  profile-fitting 
strategy.  It  can  be  simplified  by  using  more  physical  constraints.  For  example,  if  the 
temperature  non-uniformity  in  a  combustion  flow  results  mainly  from  the  non-uniform 
local  equivalence  ratio,  the  mole  fraction  distribution  of  the  absorber,  e.g.  water  vapor, 
can  be  assumed  to  be  similar  to  the  temperature  distribution,  i.e. 

&(Xcbari  ^ Xchar 9  X )  C' f^char*^Tchar*X)  (2.4.5) 

where  c  is  a  constant  to  be  inferred  from  the  fitting  along  with  Tchar  and  LTchar.  When  such 
a  relationship  is  utilized  in  the  profile  fitting,  the  number  of  unknowns  will  be  reduced 
signficantly.  As  another  example,  the  mole  fraction  can  be  assumed  to  be  constant  along 
the  measurement  path  in  cases  where  the  temperature  non-uniformity  is  much  more 
significant  than  that  of  the  mole  fraction.  Finally,  measurements  by  other  sensors,  CFD 
calculations,  past  knowledge  on  the  target  flow  field  or  a  similar  system,  etc.,  can  be  used 
to  constrain  some  of  the  variables  and  reduce  the  number  of  unknowns  in  the  postulated 
temperature  profile.  For  example,  thermocouples  can  be  used  to  measure  the  gas 
temperature  near  the  wall,  T„,  so  that  this  quantity  is  known. 

The  profile-fitting  strategy  obviously  requires  using  physical  understanding  of  the 
flow  fields  to  constrain  the  temperature  and  mole  fraction  distributions  along  the 
measurement  path.  Simulation  studies  have  demonstrated  that  the  measurement  results 
will  be  improved  significantly  by  using  as  many  physical  constraints  as  possible  to 
reduce  the  number  of  unknowns  in  the  postulated  temperature  profile. 

Principles  of  Temperature  Binning 

If  a  priori  knowledge  of  the  flow  fields  is  not  readily  accessible,  the  temperature- 
binning  method  can  be  utilized  to  extract  the  characteristic  information  on  the  LOS  non¬ 
uniformities.  The  temperature-binning  method  is  also  useful  for  cases  where  the 
temperature  distribution  of  the  flow  fields  cannot  be  represented  by  a  simple  profile  or 
cases  where  obtaining  an  exact  temperature  distribution  profile  is  not  necessary  for  the 
application.  The  temperature-binning  method  is  derived  from  a  discretization  of  the  path 
integral  of  the  absorbance 

*  =  J  =  >"WO-  (2-4.6) 

r  y=l 

which  implies  decomposing  the  entire  measurement  path  with  non-uniform  properties 
into  n  sections,  each  with  a  nearly  uniform  temperature  of  7),  absorber  mole  fraction  of 
Xabsj  and  path  length  of  Lj.  For  the  selected  m  absorption  transitions,  the  following  linear 
equation  set  can  be  established: 
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where  the  n  temperature  bins  are  prescribed  based  on  a  rough  estimation  of  the  possible 
temperature  range  along  the  measurement  path.  Once  the  number  of  absorption 
transitions  is  larger  than  the  number  of  temperature  bins,  i.e.  m>  n.  The  set  of  equations 
2.4.7  can  be  solved  by  constrained  linear  least-square  fitting 

.  (2.4.8) 


such  that  ( XabsL)j  >0  (j  =  1, •••«). 

The  solution  for  ( XabsL)j ,  called  column  density,  is  actually  the  probability  density 
function  (PDF)  of  the  absorbing  species.  More  physically  meaningful  results  can  be 
inferred  from  the  PDF  solution  if  other  physical  constraints  are  available.  For  example,  if 
the  mole  fraction  of  the  absorbing  species  can  be  assumed  to  be  constant,  the  fraction  of 
path  length  (jj)  for  each  temperature  bin  can  be  calculated  from  the  column  density  since 


f  ^ab,L)l 


(2.4.9) 


Figure  2.4.2  shows  such  a  case,  where  the  combustion  gas  temperature  along  the 
measurement  path  has  a  parabolic  distribution,  as  plotted  in  panel  a,  and  the  mole  fraction 
of  absorber  can  be  assumed  to  be  constant.  Ten  temperature  bins  within  a  typical 
combustion  temperature  range  of  1000-2000  K  have  been  prescribed.  Panel  b  shows  the 
exact  PDF  solution,  which  reveals  that  no  temperature  is  higher  than  1 800  K  along  the 
measurement  path  since  the  two  highest-temperature  bins  are  empty  and  ~  35%  of  the 
entire  path  is  at  the  temperature  of  1700-1800  K,  while  only  6%  is  at  1000-1 100  K,  etc. 
This  example  illustrates  that  non-intrusive  and  rapid  measurements  of  this  type  can  be 
useful  in  characterizing  channel  flows,  e.g.,  in  connection  with  process  control. 


T[K] 

Figure  2.4.2.  Parabolic  temperature  distribution  (panel  a)  and  the  temperature-binning  result  for  this 
distribution  (panel  b). 


34 


Simulation  studies  have  shown  that  increasing  the  number  of  bins  n  (for  fixed  number  of 
transitions  rri)  will  deteriorate  the  measurement  accuracy.  Thus,  a  moderate  number  of 
bins,  e.g.  five,  should  be  used  initially  and  more  bins  justified  only  if  more  transitions 
can  be  measured.  In  spite  of  the  limited  number  of  bins  and  the  lack  of  information  on 
the  spatial  arrangement  of  the  bins  along  the  path,  the  PDF  solution  is  sufficient  for 
many  monitoring  and  control  applications  where  the  goal  focuses  on  minimizing  or 
maximizing  the  non-uniformities,  e.g.  pattern  factor  sensing  and  control,  or  on 
monitoring  the  fluctuations  in  the  PDF  solution  in  studies  of  turbulence. 

Demonstration  Experiments 

The  two  strategies  to  interpret  LOS  absorption  measurements  in  flows  with  a  non- 
uniform  temperature  distribution  are  demonstrated  by  laboratory  measurements  of  a  “2- 
Zone”  temperature  distribution  with  a  wavelength-multiplexing  (WM)  scheme.  H2O 
vapor  is  selected  as  the  absorbing  species  because  it  is  naturally  present  in  the 
atmosphere,  it  is  a  major  product  of  hydrocarbon  combustion,  and  it  has  strong 
rovibrational  overtone  transitions  in  the  near-infrared  that  overlap  with  readily  available 
telecommunication  diode  lasers. 

A  “2-zone”  temperature  LOS  experiment  is  conducted  to  illustrate  the  measurement 
and  data  analysis  concepts.  The  schematic  of  the  experiment  is  shown  in  Fig.  2.4.3.  The 
free  space  laser  beam  has  a  total  path  length  of  54  cm  set  by  the  spacing  between  the 
pitch  lens  and  the  catch  lens.  A  25.4  cm  long  flat-flame  burner  is  located  at  the  center  of 
the  measurement  path  to  create  a  uniform  hot  section.  Therefore,  the  entire  LOS 

I  _  54.0  cm  I 


25.4  cm 
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Figure  2.43.  Schematic  of  the  experimental  setup  for  a  five  laser  multiplexed  absorption  sensor. 


measurement  path  can  be  divided  roughly  into  two  zones,  one  in  the  cold  room  air  and 
the  other  in  the  hot  flame  gases.  This  experimental  arrangement  is  illustrative  of 
combustion  flows  with  confined  flames  that  do  not  fill  the  flame  chamber  or  duct 
completely.  The  flat  flame  burner  is  fueled  by  premixed  ethylene  and  dry  air,  and  a  stable 
laminar  flame  can  be  obtained  over  the  equivalence  ratio  range  of  0.6-1 .4.  For  the  current 
demonstration  experiments,  the  ethylene  and  air  flow  rates  are  measured  to  be  1.8  slm 
(standard  liters  per  minute)  and  34.0  slm,  respectively,  by  calibrated  rotameters.  The 
equivalence  ratio  is  thus  0.76,  which  leads  to  an  equilibrium  water  vapor  mole  fraction  of 
0.1  with  an  uncertainty  of  ±0.03  due  to  the  measurement  uncertainty  of  the  fuel/air  flow 
rates. 
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The  flame  temperature  is  measured  at  the  height  of  the  laser  beam  (5  mm  above  the 
burner  surface)  by  a  type  S  thermocouple  with  a  bead  size  of  2  mil  (5 1  /mi),  and  radiation 
corrections  of  55  K  are  applied.  The  core  part  of  the  flame  (1-24  cm)  has  a  uniform 
temperature  distribution.  The  average  is  1534  K,  with  an  uncertainty  of  3%  estimated 
from  the  scatter  of  the  thermocouple  readings  and  the  uncertainty  in  the  radiation 
correction.  This  measured  flame  temperature  is  lower  than  the  adiabatic  flame 
temperature  at  the  measured  equivalence  ratio  mainly  due  to  the  radiation  loss  to  the 
surroundings  and  the  heat  conduction  to  the  water-cooled  burner  surface.  At  both  edges 
of  the  flame,  there  is  a  very  slight  temperature  rise  and  then  a  sharp  drop  to  the  room 
temperature,  thus  creating  well-defined  high-  and  low-temperature  zones  along  the  LOS 
measurement  path. 

The  thermocouple  measurement  of  the  room  temperature  is  298  K,  which  agrees  with 
the  readings  of  a  mercury  thermometer.  The  room  air  humidity  is  measured  by  a 
calibrated  hygrometer  to  be  56%,  which  indicates  a  water  mole  fraction  of  1 .75%. 

Wavelength-Multiplexed  Diode  Laser  Absorption  Measurements 

The  layout  of  the  WM  scheme  also  is  shown  in  Fig.  2.4.3.  The  fiber-coupled  outputs 
from  five  distributed-feedback  (DFB)  lasers  are  combined  and  coupled  into  one  single¬ 
mode  fiber  by  a  multiplexer  (Zolo  Technologies).  The  multiplexed  laser  beam  is 
collimated  by  the  pitch  lens  and  propagates  across  the  “2-Zone”  measurement  path. 
Another  lens  is  used  to  catch  the  free  space  laser  beam  and  focus  it  into  a  multi-mode 
fiber  that  leads  to  a  grating-based  demultiplexer  (Zolo  Technologies).  The  wavelength- 
multiplexed  beam  then  is  diffracted  into  the  constituent  five  wavelengths,  collected  on 
fiber  for  each  channel  and  delivered  to  five  InGaAs  detectors. 

The  five  DFB  lasers  emit  near  1343,  1345,  1392,  1395  and  1398  nm.  The  fiber- 
coupled,  echelle  grating  based  ZMUX™  multi/demultiplexers  are  chosen  over  free  space, 
conventionally  ruled  gratings  because  these  compact  devices  allow  dense  wavelength 
multiplexing  and  demultiplexing  (up  to  44  wavelengths)  with  high  efficiency,  pre-set 
alignment,  and  well-controlled  polarization-dependent  loss  and  thermal  drift. 

The  DFB  lasers  can  be  scanned  across  2-3  cm'1  by  injection  current  variation,  enabling 
access  to  the  seven  strong  water  vapor  transitions  listed  in  Table  1.  These  seven 
transitions  are  used  for  this  experiment  to  illustrate  the  proposed  sensor  concepts  and 
demonstrate  the  feasibility,  though  they  are  not  the  optimum  choice  for  the  current 
temperature  range.  All  five  lasers  are  scanned  simultaneously  at  I  kHz  with  a  linear 
current  ramp.  The  transmitted  laser  signals  obtained  by  the  five  detectors  are  recorded 
simultaneously  at  a  5  MHz  sampling  rate. 

Data  Reduction 

The  recorded  raw-data  scans  from  each  of  the  five  channels  are  corrected  for  detector 
DC  off-sets  and  background  emission,  although  most  of  the  emission  from  the  flame  has 
been  rejected  by  the  ZMUX™  demultiplexer,  which  acts  as  a  band-pass  filter  with  a 
narrow  bandwidth  of  ~1  nm  for  each  channel.  The  corrected  raw  data  for  each  channel 
then  are  averaged  for  every  ten  sequential  scans  to  reduce  stochastic  noise.  The  resulting 
time  response  of  this  demonstration  measurement  is  10ms,  owing  to  the  ten  sweep 
average  of  the  I  kHz  wavelength  scans.  (This  time  response  is  not  a  maximum  for  the 
technique,  as  we  have  routinely  made  scanned-wavelength  direct  absorption 
measurements  at  10kHz  and  have  performed  limited  demonstration  measurements  at 
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scan  rates  of  50  kHz  .)  From  each  averaged  laser  scan  (i.e.  the  transmitted  signal  /,),  the 
unattenuated  laser  intensity  (i.e.  the  baseline  7o)  is  determined  by  fitting  the  part  of  the  /, 
trace  without  absorption  with  a  polynomial. 

Table  1  Spectroscopic  parameters  of  the  seven  water  vapor 
transitions  used  in  the  demonstration  experiments. 


Line  Index 

vt[cm  '] 

S(T0)  [cm ^atm'1] 

1 

7154.35 

3.670E-04 

1789.04 

2 

7164.90 

3.550E-03 

1394.81 

3 

7185.60 

1 .960E-02 

1045.06 

4 

7417.82 

1.070E-02 

1079.08 

5 

7444.36 

1 . 1 00E-03 

1786.00 

6 

7164.07 

2.021E-02 

300.36 

7 

7419.17 

2.222E-02 

842.36 

The  integrated  absorbances  A,  for  the  seven  lines  are  calculated  from  the  measured 
absorption  spectra.  The  data  reduction  is  more  complicated  than  that  for  uniform  cases 
since  the  lineshape  measured  along  a  non-uniform  temperature  distribution  no  longer  can 
be  modeled  by  a  single  Voigt  function.  Once  the  integrated  absorbances  for  the  seven 
lines  are  obtained,  either  the  profile-fitting  strategy  or  the  temperature-binning  strategy 
can  be  applied  to  characterize  the  non-uniform  temperature  distribution  using  the  models 
above.  The  spectroscopic  parameters  of  the  seven  transitions  used  in  the  calculations  are 
listed  in  Table  1. 

If  we  assume  that  both  the  temperature  and  water  mole  fraction  along  the  LOS 
measurement  path  are  uniform,  this  “uniform”  temperature  can  be  inferred  from  the  ratio 
of  measured  absorbance  of  any  line  pair.  Here  line  pair  I  and  2  yields  1 199  K,  while  line 
pair  1  and  3  produces  1122  K.  This  discrepancy  suggests  non-negligible  non- 
uniformities,  thus  illustrating  the  value  of  the  new  strategy  to  use  multiple  absorption 
measurements  to  characterize  the  actual  temperature  distribution. 

Profile-Fittine  Results 

In  the  profile-fitting  calculation,  the  shape  of  the  non-uniform  property  distribution 
must  be  postulated  in  advance.  The  layout  of  the  experimental  setup  (hot  flame 
temperature  in  the  middle  with  cold  room  temperatures  on  both  sides)  enables  us  to 
postulate  a  “2-Zone”  profile,  as  shown  in  Fig.  2.4.4,  to  model  the  temperature  and  mole 
fraction  distribution  along  the  LOS  measurement  path.  The  unknowns  (free  parameters) 
to  be  solved  from  the  multiple  absorption  measurements  are  the  temperature  and  mole 
fraction  for  each  zone,  i.e.  T\,  T2,  X\  and  X2.  We  investigate  four  different  interpretations 
of  the  measured  absorption  data,  each  with  a  different  number  of  constraints  on  the 
postulated  distributions.  For  each  of  these  four  cases,  the  influence  of  the  number  of  lines 
on  the  sensor  performance  is  investigated. 


37 


Figure  2.4.4.  The  “2-Zone”  property  distribution  postulated  for  profile-fitting  calculation. 

In  case  1,  the  measured  absorbances  are  fit  with  all  four  unknowns  T\,  T2,  X\  and  X2 
as  free  parameters,  and  the  average  values  are  listed  in  Table  2.  Both  the  measurement 
accuracy,  as  indicated  by  the  difference  between  the  measured  and  the  expected  values, 


Table  2  The  average  values  of  the  profile-fitting  results  (and 
percentage  deviation  from  the  expected  values)  with 
different  number  of  lines  for  cases  1  and  2. 


Lines 

Case  1 

Case  2 

T.PC] 

7-2[K] 

*.[%] 

X2[%] 

7-,[K] 

7-2  [K] 

Xi[%] 

309 

1374 

1.1 

8.4 

227 

1264 

7.8 

1-5 

(4%) 

(-10%) 

(-37%) 

(-16%) 

(-24%) 

(-18%) 

(-22%) 

292 

1489 

1.7 

9.4 

295 

1524 

9.7 

1-7 

(-2%) 

(-3%) 

(-3%) 

(-6%) 

(-1%) 

h%) 

(-3%) 

Expected 

298 

1534 

1.75 

10.0 

298 

1534 

10.0 

and  the  statistical  precision,  increase  appreciably  with  the  addition  of  data  from  lines  6 
and  7.  By  using  the  absorbances  measured  on  lines  1-5,  the  temperature  and  mole- 
fraction  inferred  deviate  significantly  from  the  expected  values,  which  is  due  partly  to  the 
narrow  span  of  the  lower  state  energy  E”  of  these  5  transitions.  By  adding  measurements 
from  line  6  and  7  with  much  lower  E”,  the  results  for  7i,  T2,  X\  and  X2  quickly  converge 
to  the  expected  values  within  an  accuracy  of  2%,  3%,  3%  and  6%  respectively.  This 
result  illustrates  the  importance  of  optimized  selection  of  absorption  transitions. 

In  case  2,  the  H2O  mole  faction  of  the  cold  zone  X\  is  fixed  at  1.75%,  which  is 
determined  using  the  measured  humidity  and  an  estimated  temperature  (300  K)  of  the 
room  air.  The  three  remaining  unknowns  T\,  T2,  and  X2  then  are  allowed  to  vary  in  the 
nonlinear  least  square  fit,  and  the  average  values  are  listed  in  Table  2.  Again  the 
measurement  accuracy  increases  systematically  with  the  addition  of  data  from  lines  6  and 
7.  By  using  all  7  transitions,  the  T\,  T2,  and  X2  can  be  measured  within  an  accuracy  of  I  %, 
1%,  and  3%,  respectively,  which  improves  from  the  corresponding  results  for  case  1. 
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Generally  speaking,  the  addition  of  physical  constraint(s)  improves  the  accuracies  of  the 
fit  values. 

In  case  3,  the  H20  mole  fractions  for  the  two  zones  X,  and  X2  are  fixed  at  1 .75%  and 
10%  (the  equilibrium  value  for  the  measured  equivalence  ratio),  respectively.  Since  there 
are  only  two  remaining  unknowns,  T\  and  T2,  a  minimum  number  of  3  lines  is  required 
for  the  nonlinear  least  squares  fit,  and  the  average  values  are  listed  in  Table  3.  These  fit 
results  reveal  the  same  trend,  namely  an  increasing  of  the  accuracy  of  the  fit  values  with 
an  increased  number  of  lines.  Even  with  the  minimum  number  of  lines,  however, 
satisfactory  results  of  7i  and  T2  can  be  obtained  (with  an  accuracy  of  4%  and  1%, 
respectively)  by  fixing  the  H2Q  mole  fraction  values. 


Table  3  The  average  values  of  the  profile-fitting  results  (and 
percentage  deviation  from  the  expected  values)  with  different 
number  of  lines  for  cases  3  and  4. 


Lines 

Case  3 

Case  4 

7i[  K] 

7-2[K] 

Ti  [K] 

Xi  [%] 

1-3 

286  (-4%) 

1516  (-1%) 

1487  (-3%) 

9.4  (-6%) 

1-5 

302(1%) 

1554(1%) 

1520  (-1%) 

9.6  (-4%) 

1-7 

299  (0%) 

1552(1%) 

1537(0%) 

9.8  (-2%) 

Expected 

298 

1534 

1534 

10.0 

In  case  4,  the  properties  of  the  room  air  T,  and  X,  are  fixed  in  the  fit,  and  thus  only  the 
properties  of  the  hot  flame  zone  T2  and  X2  are  the  free  parameters,  and  the  average  values 
are  listed  in  Table  3.  The  accuracy  increases  with  the  number  of  lines,  particularly  with 
the  inclusion  of  lines  6  and  7.  Using  all  7  transitions,  T2  and  A^are  determined  within  an 
accuracy  of  1  %  and  2%,  respectively. 

This  two-zone  measurement  arrangement  clearly  is  applicable  to  cases  with 
interference  absorption  by  humid  room  air  bounding  a  relatively  uniform  hot  zone  of 
interest.  Only  the  temperature  and  water  content  of  the  hot  target  gas  typically  are 
desired,  but  sometimes  room  air  inevitably  is  enclosed  in  the  LOS  beam  path  near  the 
pitch  or  catch  optics,  since  the  laser  sources  or  the  measurement  geometries  prohibit 
using  fibers  to  deliver  the  laser  beam  directly  to  the  edge  of  the  flow  fields  to  be 
measured.  Normally,  in  our  laboratory,  the  open  path  in  this  boundary  is  purged  with  N2 
to  remove  the  interference  absorption  by  the  room  air,  but  the  purging  efficiency  cannot 
be  guaranteed  for  all  applications,  especially  if  the  purging  system  is  not  well-designed, 
installed,  and  maintained.  By  using  more  than  two  transitions,  as  is  done  in  case  4,  the 
flow  field  properties  of  interest  can  be  inferred  directly  from  the  LOS  absorption  data 
without  installing  complex  purging  systems. 

Table  4  shows  a  comparison  of  the  profile-fitting  results  for  all  four  cases  when  data 
from  all  7  transitions  are  used.  It  confirms  that  the  fitting  results  improve  significantly 
when  the  number  of  unknowns  (free  parameters)  in  the  postulated  distribution  profiles 
can  be  reduced  by  using  physical  constraints. 

Temperature-Binning  Results 

Alternatively,  we  can  interpret  the  absorption  data  using  the  temperature-binning 
strategy.  An  estimation  of  the  possible  temperature  range  along  the  LOS  measurement 
path  allows  the  prescription  of  five  temperature  bins  as  shown  in  Fig.  2.4.5.  All  seven 
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Table  4  Comparison  of  the  profile-fitting  results  (and  percentage 


deviation  from  the  expected  values)  for  all  four  cases  with 
all  7  lines. 


Case 

r.pq 

7HK] 

X\  [%] 

x2[%] 

1 

292  (-2%) 

1489  (-3%) 

1.7  (-3%) 

9.4  (-6%) 

2 

295  (-1%) 

1524  (-1%) 

~ 

9.7  (-3%) 

3 

299  (0%) 

1552(1%) 

— 

— 

4 

— 

1537  (0%) 

— 

9.8  (-2%) 

Expected 

298 

1534 

1.75 

10.0 

transitions  are  used  in  the  analysis.  The  resultant  column  densities  (Xmo L)j,  as 
represented  by  the  solid  bars,  are  very  close  to  the  expected  PDF  solutions,  as  indicated 
by  the  dash-dot  lines.  Only  the  two  side  bins  have  non-zero  solutions.  They  suggest  that 
the  LOS  measurement  path  can  be  modeled  approximately  as  two  zones,  one  at  -300  K 
and  the  other  at  -1500  K,  which  are  very  close  to  the  conditions  of  the  experiment.  The 
population  of  water  vapor  in  the  highest  temperature  bin  is  much  larger  than  that  in  the 
lowest  temperature  bin,  which  implies  that  the  “2-T”  distribution  might  result  from  a  cold 
spot  somewhere  along  an  otherwise  uniform  high  temperature  region  or  a  cold  boundary 
layer  on  both  sides  of  a  uniform  hot  core  temperature.  Up  to  this  point,  no  a  priori 
knowledge  of  the  flow  fields,  except  for  an  estimate  of  the  possible  temperature  range, 
has  been  applied  to  extract  the  characteristic  information  on  the  LOS  non-uniformities 
discussed  above. 
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Figure  2.4.5.  Illustration  of  the  temperature-binning  results  solved  using  all  7  transitions  with  five  bins. 


Current  Status 

We  have  performed  an  initial  laboratory  demonstration  of  two  strategies  for 
measurements  of  non-uniform  temperature  distributions  using  Iine-of-sight  (LOS) 
absorption  spectroscopy.  Based  on  the  measurements  of  multiple  absorption  transitions 
with  different  temperature  dependences,  the  profile-fitting  strategy  fits  a  postulated 
temperature  distribution  profile,  while  the  temperature-binning  strategy  determines  the 
temperature  probability  density  function  (PDF)  along  the  LOS  using  prescribed 
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temperature  bins.  Both  strategies  can  be  modeled  mathematically  as  least-square  fitting 
problems.  In  the  non-optimized  demonstration  experiment,  a  “2-Zone”  temperature/mole- 
fraction  distribution  created  in  laboratory  is  measured  with  a  wavelength-multiplexing 
scheme  where  five  DFB  diode  lasers  are  multiplexed  to  probe  seven  absorption 
transitions  of  water  vapor.  The  absorption  data  are  analyzed  by  both  strategies  to 
illustrate  the  sensor  concepts  and  investigate  the  sensor  performance.  The  experimental 
results  demonstrate  that  a  non-uniform  temperature  distribution  can  be  characterized  with 
either  strategy  by  measuring  the  LOS  absorption  for  a  limited  number  of  transitions  with 
different  temperature  dependences.  The  measurement  accuracies  will  increase  with  the 
number  of  transitions  and  with  the  use  of  optimally  chosen  transitions.  Use  of  known 
physical  constraints  also  provides  important  value  in  properly  interpreting  measurement 
results  and  improving  the  sensor  performance.  Further  work  is  planned  to  develop  design 
rules  for  the  selection  of  absorption  transitions.  Using  an  optimized  set  of  absorption 
transitions,  the  next  step  envisioned  includes  demonstration  experiments  in  systems  with 
realistic  temperature  profiles. 
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2.4.2.  Scanned-wavelength-modulation  for  high  SNR  absorption 

A  wavelength-multiplexed  diode  laser  absorption  sensor  is  developed  and  tested  for 
the  measurement  of  gas  temperature  and  the  concentration  of  water  vapor  in  a  model 
scramjet  combustor  in  Test  Cell  22  of  the  Air  Force  Research  Laboratory  (AFRL) 
Propulsion  Directorate  at  Wright-Patterson  AFB.  This  sensor  uses  multiple  diode  lasers 
to  monitor  beam  steering  and  three  spectral  features  of  water  vapor  in  the  near  infrared 
region.  The  measurements  are  made  near  the  exit  of  the  scramjet  combustor.  An 
important  aspect  of  the  work  is  the  optical  engineering  used  to  overcome  beam  steering 
from  the  high-speed  flow  and  fiber-mode-noise.  Test  measurements  at  AFRL  are  made 
with  the  sensor  using  two  detection  strategies:  (1)  scanned-wavelength  direct  absorption 
and  scanned-wavelength-modulation  spectroscopy  (scanned-WMS)  with  second 
harmonic  detection  (2 \f),  and  the  noise  performance  of  the  two  sensor  architectures  are 
compared.  The  WMS-2f  technique  produces  superior  noise  rejection  and  improved 
overall  signal  to  noise  ratios  compared  to  direct  absorption. 

Background 

Traditionally,  scanned-wavelength  direct  absorption  spectroscopy  has  been  the 
technique  of  choice  for  aero-engine  applications  because  of  its  simplicity,  accuracy,  and 
ability  to  make  absolute  measurements  [2.4.1],  In  harsh  combustion  environments  such 
as  the  scramjet,  however,  the  technique  can  suffer  from  low  signal-to-noise  ratios  arising 
from  mechanical  and  gasdynamic  beam  steering,  weak  absorption  at  high  temperatures, 
and  fiber  mode  noise.  The  scanned-direct  absorption  technique  must  rely  entirely  on 
careful  optical  engineering  to  reduce  noise  to  acceptable  levels.  Another  technique  that 
has  been  explored  for  its  increased  sensitivity  to  small  signals  and  ability  to  reject  noise  is 
wavelength  modulation  spectroscopy  (WMS)  with  second  harmonic  detection  (2 \f).  This 
technique,  often  simply  referred  to  as  WMS-2/  employs  a  lock-in  amplifier  to  isolate  the 
second  harmonic  of  the  detected  signal,  greatly  enhancing  the  noise  rejection 
characteristics.  Example  signal-to-noise  ratio  (SNR),  temperature,  and  fast-Fourier 
transform  (FFT)  data  obtained  in  a  scramjet  combustor  using  both  scanned-direct 
absorption  and  WMS-2/ under  nearly  identical  operating  conditions  will  be  compared  to 
illustrate  the  benefits  of  WMS-2/in  harsh  environments. 

Sensor  Architecture  and  Experimental  Design 

The  sensor  utilizes  a  beam  of  infrared  light  consisting  of  the  multiplexed  output  of 
four  TDLs.  The  light  is  pitched  across  the  combustor  exit  plane,  collected,  and  dispersed 
onto  four  separate  detectors.  Three  of  the  lasers  are  wavelength-scanned  across  H2O 
transitions  and  used  to  infer  the  static  temperature  via  direct  absorption  and  WMS-2f. 
The  fourth  laser  is  fixed  at  a  non-resonant  wavelength  to  detect  transmission 
perturbations  from  gasdynamic  or  mechanical  beam  steering,  window  fouling,  or 
particulates  in  the  flow. 

A  schematic  of  the  optical  arrangement  of  the  sensor  is  shown  in  Fig.  2.4.6.  Four 
fiber-coupled  DFB  diode  lasers  are  multiplexed  with  a  9-micron  single-mode  fiber 
combiner.  A  10-meter  single-mode  fiber  delivers  the  beam  to  the  scramjet  test  section, 
where  it  is  pitched  across  the  combustor  exit  with  an  aspheric  collimating  lens.  After 
traversing  the  test  region,  the  beam  is  collected  with  a  large  diameter  lens  and  focused 
into  a  400  micron  multimode  fiber.  The  beam  travels  through  the  multimode  fiber  to  a 
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demultiplexing  setup  a  few  meters  from  the  test  section.  The  light  exits  the  multimode 
fiber  and  is  collimated  by  a  large  diameter  aspheric  lens  onto  a  30  mm  square,  1200 
grooves/mm  diffraction  grating.  The  dispersed  light  beams  are  focused  onto  separate 
3mm  diameter  InGaAs  photodetectors. 
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Figure  2.4.6.  Schematic  of  optical  setup. 


The  scramjet  test  facility  at  Wright-Patterson  Air  Force  Base  in  Dayton,  Ohio  will  not 
be  covered  in  detail  here  because  an  in-depth  treatment  is  available  in  the 
literature[2.4.3].  Compressed  air  (up  to  13.6  kg/s  and  920K)  is  vitiated  to  raise  the  pre¬ 
expansion  temperature.  Make-up  oxygen  replaces  that  burned  in  the  vitiator  heater.  The 
gas  stream  is  expanded  to  supersonic  speeds  through  a  Mach  2.7  nozzle.  Liquid  JP-7  fuel 
is  heated  to  a  vaporized  condition  and  injected  to  the  near  atmospheric  pressure  gases  just 
upstream  of  a  cavity-based  flameholder. 

The  combustor  is  fitted  with  2.5  cm  x  12  cm  quartz  optical  access  ports  located  16.5 
cm  upstream  of  the  end  of  the  combustor.  The  outward  face  of  each  window  is  wedged 
to  avoid  optical  interference  effects  that  arise  when  scanning  laser  wavelength  in  the 
presence  of  parallel  optical  faces  in  the  beam  path.  The  collimating  pitch  and  catch 
optics  are  arranged  on  translation  stages  that  can  traverse  the  optical  ports  vertically 
while  maintaining  beam  alignment. 

All  absorption  spectroscopy  techniques  rely  on  careful  selection  of  the  quantum 
transitions  (lines)  that  will  be  used  to  probe  the  gaseous  flow.  Careful  selection  begins  by 
choosing  transitions  of  a  species  present  in  the  flow  (in  this  case,  water  vapor)  with 
desirable  absorption  properties  for  the  particular  application.  For  these  experiments, 
three  lines  are  chosen  with  different  lower  state  energies.  The  lower  state  energy, 
determines  the  equilibrium  molecular  population  in  the  absorbing  state  as  a  function  of 
temperature  and  thus  can  be  used  to  describe  how  the  absorption  strength  of  a  particular 
transition  varies  with  temperature.  The  three  transitions  used  for  these  experiments  are 
shown  in  Table  2.4.5.  The  lower  state  energies  are  chosen  to  give  peak  temperature 
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sensitivity  for  the  1392  nm  and  1343  nm  features  at  500  K  and  the  1392  and  1469  nm 
features  at  1600  K.  These  temperatures  coincide  with  the  expected  temperature  of  a 
scramjet  run  with  only  the  vitiator  operating  and  with  full  combustion,  respectively. 

The  three  lines  also  are  chosen  for  their  isolation  from  other  transitions,  laser 
availability,  and  correct  wavelength  spacing  to  allow  dispersal  by  the  grating.  The  key 
spectral  parameter  of  line  strength  is  validated  for  each  feature  under  carefully  controlled 
conditions  in  a  laboratory  furnace. 


Table  2.4.5.  Line  strength  and  lower  state  energy  parameters  for  spectroscopic 
_ _ features  used  in  these  experiments. _ 


Wavelength 

nm 

(Hitran  ’04) 

Frequency 

cm'1 

(Hitran  '04) 

Line  Strength  (296  K) 
cm'2  atm'1 
(Hitran  '04) 

Line  Strength  (296  K) 
cm'2  atm'1 
(Measured) 

Lower  State  Energy 
cm'1 

(Hitran  *04) 

1392 

7185.6 

1.97(10)"z 

1.96(10)"z 

1045.1 

1343 

7444.35  and  7444.37 

1.12(H))-3 

I.IO(IO)-3 

1774.8  and  1806.7 

1469 

6807.83 

1.02(1 0)^ 

6.48(10)-' 

3319.4 

Measurement  Techniques:  Scanned-Wavelength  Direct  Absorption 


Scanned-wavelength  direct  absorption  has  been  a  staple  of  TDL  sensors  for  harsh 
environments.  The  important  aspects  of  the  technique  are  shown  in  Fig.  2.4.7.  The 
injection  current  of  the  TDL  is  driven  with  a  linear  ramp  signal  (in  this  case,  a  4  kHz 
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Figure  2.4.7.  (top)  Raw  detector  signal  during  a  full  combustion  run  for  a  single  direct  absorption 
laser  scan  with  baseline  fit  (1392  nm  feature);  (bottom)  Absorbance  plot  resulting  from  laser  scan  and 
baseline. 

sawtooth  signal).  The  linear  ramp  in  injection  current  produces  a  near-linear  variation  in 
laser  intensity  that  can  be  seen  in  the  top  panel.  The  variation  in  laser  intensity  is 
accompanied  by  a  near-linear  change  in  laser  wavelength.  As  the  laser  is  tuned  across  the 
absorption  feature  some  light  is  absorbed,  reducing  the  detector  signal.  A  baseline  fit  is 
performed  on  the  non-absorbing  “wings”  of  the  feature  to  infer  the  laser  intensity  in  the 
absence  of  absorption.  The  Beer-Lambert  relation  (Eqn.  2.4.10)  can  be  rearranged  to 
yield  Eqn.  2.4.1 1,  which  gives  the  absorbance  (a)  as  a  function  of  the  transmitted  laser 
intensity  (It,  the  detector  signal),  and  the  incident  laser  intensity  (I,*  the  baseline  fit).  The 
absorbance  plot  resulting  from  the  laser  scan  is  shown  in  the  lower  panel. 
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(2.4.10) 


—  =  exp(-a) 

K 


(2.4.11) 


The  absorbance  is  a  function  of  line  strength  (S),  species  partial  pressure  (Pj), 
pathlength  (L),  and  line  shape  function  (<j>),  as  shown  in  Eqn.  2.4.12. 


a  -  S(T ) xPtxLx <f>(v) 


(2.4.12) 


The  line  shape  is  a  normalized  function  that  describes  the  profile  of  a  spectroscopic 
transition  in  frequency  space.  Thus  if  the  measured  absorbance  of  a  transition  is  fit  with 
the  proper  profile  (a  Voigt  profile  is  used  here)  and  integrated  with  respect  to  frequency, 
the  line  shape  function  is  equal  to  unity  and  is  no  longer  important.  This  step  is  shown  in 
Eqn.  2.4.13,  where  the  result  from  integrating  with  respect  to  frequency  is  called  the 
integrated  absorbance  area. 

00 

area  -  S(T) x  Pfx  Lx  J^(v) •  dv  =S(T) xP(xL  (2.4.13) 


For  thermometry,  the  ratio  of  the  integrated  area  of  two  absorption  features  is  used. 
When  the  ratio  is  calculated,  the  species  partial  pressure  and  pathlength  cancel.  The 
remaining  ratio  of  line  strengths  is  purely  a  function  of  temperature  and  is  used  to  infer 
the  temperature  from  the  measured  area  ratio. 


area 


area 


S\(T) 


2  S2(T ) 


(2.4.14) 


The  species  concentration  is  calculated  using  Eqn.  2.4.13  with  one  spectral  feature.  The 
area  and  pathlength  are  known,  and  the  line  strength  is  calculated  using  the  temperature 
previously  inferred  from  the  ratio. 


Measurement  Techniques;  Scanned- WMS-2f 

The  theory  of  WMS  with  2/  detection  has  been  covered  extensively  in  the  literature 
[2.4.4],  It  is  used  generally  in  situations  involving  small  absorbances  and/or  noisy 
environments.  In  its  most  basic  form,  WMS-2/  involves  rapidly  modulating  the  laser 
wavelength,  passing  the  light  through  an  absorbing  medium,  and  running  the  detector 
signal  through  a  lock-in  amplifier  set  to  isolate  the  signal  at  twice  the  modulation 
frequency.  The  resulting  signal  can  be  used  to  calculate  various  gasdynamic  properties. 
In  the  scanned-wavelength  variant  of  WMS-2/  employed  here,  the  rapidly  modulated 
laser  wavelength  also  is  tuned  slowly  across  an  absorption  feature,  and  the  peak  height  of 
the  WMS-2/ signal  is  used  to  infer  gas  properties.  The  top  panel  of  Fig.  2.4.8  depicts  the 
raw  detector  signal  of  a  single  sweep  of  scanned- WMS-2/  As  with  the  direct  absorption 
method  discussed  before,  the  laser  injection  current  is  driven  with  a  linear  ramp  (2  kHz), 
but  now  there  is  also  a  high  frequency  sinusoid  superimposed  on  the  ramp  (190  kHz). 
The  laser  intensity  varies  with  the  drive  signal  and  is  accompanied  by  a  similar  variation 
in  laser  wavelength. 
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Figure  2.4.8.  (top)  Raw  detector  signal  during  a  full  combustion  run  for  a  single  laser  scan  with 
modulation  (1392  nm  feature).  Also  shown  is  fit  to  obtain  G*I0(t).  (bottom)  WMS-2/  signal  resulting 
from  laser  scan. 


The  laser  frequency  (inverse  wavelength),  v(t),  and  the  laser  intensity,  I0(t),  for  the 
situation  described  above  are  given  by  Eqn.  2.4.15  and  2.4.16. 

v(/)  =  v(/)  +  a  cos(cot)  1 5) 

hi  0  =  4(0D  +  i0(t)cos(cot  +  if/)]  (2.4.16) 


where  v(/)and  lQ(t)  are  the  average  laser  frequency  and  intensity  that  are  varying  slowly 

due  to  the  linear  ramp,  co  is  the  modulation  frequency,  a  is  the  frequency  modulation 
amplitude,  and  /<,(/)  is  the  intensity  modulation  amplitude  (that  is  a  constant  normalized 
by  /.(#)  )•  V  is  the  phase  shift  between  frequency  and  intensity  and  is  commonly 
assumed  in  the  literature  to  be  71  radians.  Due  to  the  small  modulation  depths  utilized 
for  this  work,  this  assumption  is  made  with  minimal  resulting  error;  however  real  diode 
lasers  have  a  value  of  <//  that  is  somewhat  greater  than  7t  and  can  be  measured  to 
investigate  the  magnitude  of  this  simplification. 

When  modulated  laser  light  at  a  frequency  corresponding  to  a  quantum  transition  is 
passed  through  an  absorbing  medium,  it  is  attenuated.  The  intensity  of  light  that  remains 
after  having  passed  through  the  medium  (/,)  is  proportional  to  a  transmission  coefficient 
(r)  that  is  a  function  of  the  frequency  of  the  light: 

/,(/)  =  I0(t)  r(v  +  a  cos  cot)  (2.4.17) 

Because  the  transmission  coefficient  is  a  periodic,  even  function  in  cot,  it  can  be  expanded 
in  a  Fourier  cosine  series.  The  mathematical  details  will  not  be  shown  here,  but  the 
resulting  signal,  normalized  by  the  average  laser  intensity  and  detector  gain  ( GIa ),  is 
given  by  Eqn.  2.4.18.  This  equation  assumes  linear  frequency  and  intensity  modulation 
with  a  phase  shift  of  n.  Equation  2.4.18  could  be  written  as  a  function  of  time  because 
the  average  laser  frequency  varies  in  time  with  the  linear  ramp  of  injection  current.  It 
will  be  advantageous  to  express  the  WMS-2/ signal  in  terms  of  average  laser  frequency 
for  the  purpose  of  comparing  with  simulations. 
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(2.4.18) 


When  taking  the  ratio  of  WMS-2/  peak  height  signals  from  two  spectral  lines  the 
partial  pressure  and  pathlength  cancel,  leaving  only  the  line  shape  integral  and  line 
strength  for  each  line.  Fortunately,  the  modulation  depth  of  each  laser  can  be  chosen 
such  that  the  line  shape  integral  is  nearly  constant  for  moderate  changes  in  pressure  and 
mixture  composition.  This  important  point  is  discussed  in  detail  by  Liu  et  al.  [2.4.5]. 
Thus,  if  WMS-2/ peak  height  simulations  are  performed  at  the  nominal  expected  partial 
pressure  for  both  lines  as  a  function  of  temperature,  the  resulting  ratio  of  the  two 
simulations  will  form  a  map  between  WMS-2/ peak  ratio  and  temperature.  This  map  will 
be  unaffected  by  moderate  changes  of  mixture  composition  and  pressure  about  the 
nominal  value  used  for  simulation. 

To  infer  gas  properties,  the  measured  signals  are  compared  with  simulations.  For 
temperature,  the  measured  WMS-2/  peak  ratio  is  compared  with  the  map  described 
above.  To  calculate  partial  pressure  of  the  species,  the  simulation  for  one  spectral  line  (at 
the  temperature  inferred  by  the  two  line  ratio)  is  ratioed  with  the  measured  peak  height  of 
that  line.  The  WMS-2/  signal  is  directly  proportional  to  the  partial  pressure  of  the 
species,  so  the  ratio  of  the  measured  WMS-2/ signal  and  the  simulation  will  exactly  equal 
the  ratio  of  the  measured  partial  pressure  and  the  partial  pressure  used  in  simulation. 


yGI°(vtx<*))measured 
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(2.4.19) 


The  method  above  yields  absolute  measurements  of  temperature  and  concentration  so 
long  as  rough  values  of  pressure  and  mole  fraction  are  known.  Thus,  careful 
normalization  of  the  WMS-2/ peak  height  and  the  use  of  full  simulations  with  judiciously 
chosen  and  measured  laser  parameters  (a,  i'0),enables  quantitative  measurements  without 
the  traditional  WMS-2/caIibration. 

The  scanned-wavelength  direct  absorption  method  works  well  in  environments  where 
spectrally-isolated  features  with  large  absorbances  can  be  used  and  noise  eliminated  with 
optical  engineering.  However,  as  is  the  case  with  scramjets  and  other  applications  to 
harsh  flows,  some  level  of  noise  always  will  be  present,  and  high  temperatures  will 
require  the  use  of  weakly  absorbing  features.  WMS-2/offers  certain  advantages  that  can 
address  these  issues. 

Baseline  fitting  is  a  key  component  of  scanned-wavelength  direct  absorption  data 
reduction;  however,  several  factors  increase  the  difficulty  of  making  a  good  Fit.  Fiber 
mode  noise  causes  time  and  wavelength-dependent  fluctuations  that  bend  and  distort  the 
intensity  ramp  of  the  detector  signal.  Gasdynamic  and  mechanical  beam  steering  causes 
signal  fluctuations,  resulting  in  further  errors.  In  the  limit  that  optical  engineering  can 
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control  the  aforementioned  noise  sources,  a  difficulty  that  is  still  present  is  interference 
from  neighboring  spectral  features.  Determining  where  baseline  fits  should  begin  and 
end  to  avoid  the  wings  of  the  main  and  neighboring  spectral  features  greatly  complicates 
baseline  fitting.  This  effect  is  apparent  in  Fig.  2.4.7  and  especially  is  evident  for  the  1343 
and  1469  nm  features  in  Fig.  2.4.9. 

For  small  modulation  depths,  the  WMS-2/ signal  in  the  absence  of  absorption  is  zero 
regardless  of  laser  intensity.  The  WMS-2/  signal  thus  rides  on  a  zero  background  and 
does  not  depend  on  any  baseline  fitting.  The  measured  WMS-2/  peak  height  must  be 
normalized  by  Gl„  in  order  to  compare  with  the  simulations.  A  low-order  polynomial  is 
fit  through  the  high  frequency  modulation  to  perform  the  normalization,  as  shown  in  the 
top  panel  of  Fig.  2.4.8.  While  it  may  seem  that  an  error  in  the  fit  through  the  WMS-2/ 
signal  would  induce  the  same  problems  as  an  error  in  the  baseline  Fit  of  a  direct 
absorption  scan,  there  is  a  difference  of  relative  measurement  scales  that  reduces  the 
relative  contribution  to  the  signal  uncertainty.  As  an  example,  a  1%  error  in  the  fit  of  Gl„ 
used  to  normalize  the  WMS-2/ signal  will  result  in  only  a  1%  error  in  the  normalized 
WMS-2/  peak  height.  However,  since  the  absorbance  is  measured  relative  to  the 
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Figure  2.4.9.  Comparison  of  direct  absorption  and  WMS-2f  signals  during  a  full  combustion  run. 

baseline  fit  of  a  direct  absorption  scan,  a  1%  error  on  the  baseline  fit  for  a  transition  with 
10%  peak  absorbance  will  cause  a  10%  error;  whereas  a  1%  error  in  the  estimate  of  lo  for 
WMS-2/ produces  a  1%  error  in  the  signal.  For  this  reason  WMS-2/ is  less  prone  to 
baseline  fitting  errors  than  direct  absorption. 

With  direct  absorption  spectroscopy  one  must  rely  entirely  on  optical  engineering  to 
reduce  noise.  WMS-2/  however,  employs  the  use  of  a  lock-in  amplifier  that  serves  to 
shift  the  detection  bandwidth  above  most  noise  sources.  The  lock-in  amplifier  multiplies 
the  detector  signal  by  a  sinusoidal  reference  wave  at  twice  the  modulation  frequency. 
The  components  of  the  detector  signal  thus  are  shifted  in  frequency  space  such  that  the 
signal  at  twice  the  modulation  frequency  becomes  a  DC  signal.  This  signal  then  is 
filtered  low-pass  to  eliminate  the  shifted  signals  outside  the  low  pass  filter  cutoff,  thereby 
eliminating  nearly  all  but  the  signal  at  twice  the  modulation  frequency.  For  the 
experiments  presented  here,  the  modulation  frequency  is  190  kHz.  The  lock-in  amplifier 
multiplies  the  detector  signal  by  a  380  kHz  sine  wave  and  then  low-pass  filters  the  result 
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with  a  75  kHz  cutoff,  8th  order  HR  cascade  filter.  The  resulting  signal  thus  only  contains 
signals  that  are  originally  at  380  kHz  +/-  75  kHz.  Because  most  noise  sources  are  below 
300  kHz  (including  1/laser  noise),  they  are  filtered  out  and  do  not  contribute  to  the  final 
signal. 

Comparison  of  Direct  Absorption  and  WMS-2f  Results 

Results  of  gas  temperature  and  H2O  concentration  are  obtained  in  the  scramjet  on 
consecutive  runs  with  the  same  combustor  conditions  to  compare  direct  absorption  and 
WMS-2/  Figure  2.4.9  shows  an  example  scan  for  each  line  obtained  during  a  full 
combustion  run  with  both  methods.  The  top  panel  shows  absorbance  data  that  come  from 
fitting  a  baseline  to  the  raw  laser  scan.  The  neighboring  features  reduce  the  percentage  of 
each  scan  that  can  be  fit  with  the  critical  zero  absorption  baseline.  Also,  the  features 
overlap,  complicating  Voigt  fitting  and  the  determination  of  the  integrated  absorbance 
area.  The  bottom  panel  shows  the  WMS-2/  data  under  the  same  scramjet  conditions. 
The  neighboring  features  appear  as  perturbations  next  to  the  main  feature  (that  are  not 
factored  into  the  SNR).  These  features  do  not  affect  the  measurement  of  the  peak  height 
because  the  zero  baseline  of  the  WMS-2/  signal  is  inherent  and  does  not  arise  from  a 
baseline  fit.  The  contribution  of  the  neighboring  features  to  the  WMS-2/ peak  height  is 
small  and  accounted  for  by  the  simulations.  Comparison  of  the  SNR  for  each  feature 
shows  the  improvements  possible  with  WMS-2/  The  SNR  for  the  high  temperature  1469 
nm  feature  is  improved  by  a  factor  of  two. 

Figure  2.4.10  shows  a  comparison  of  the  measured  temperatures  using  the  1392/1469 
nm  line  pair  with  direct  absorption  and  WMS-2/  The  WMS-2/  data  have  a  2  kHz 
sampling  rate  (1  WMS-2/ peak  per  laser  scan  at  2  kHz),  so  the  direct  absorption  data  are 
filtered  to  2  kHz  from  their  original  4  kHz  sampling  rate  to  give  a  fair  comparison.  The 
figure  shows  an  improvement  in  the  standard  deviation  of  the  temperature  fluctuations 
from  90  K  to  50  K. 
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Figure  2.4.10.  Comparison  of  measured  temperature  using  direct  absorption  and  WMS-2/  with  the 
1392/1469  nm  line  pair.  The  direct  absorption  data  is  filtered  from  4  kHz  to  2  kHz  to  match  the  WMS- 
2f  sampling  rate.  <|>  =  1 .02. 


Figure  2.4.11  shows  0.75  seconds  of  representative  temperature  and  H2O  partial 
pressure  data.  A  portion  of  the  start-up  transient  of  the  scramjet  combustor  is  shown. 
The  mean  temperature  increases  for  approximately  0.45  seconds  before  stabilizing,  while 
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the  mean  H2O  partial  pressure  remains  stable  throughout.  These  data  can  be  used  to 
characterize  transients  in  the  scramjet  that  may  be  difficult  to  predict  with  models. 
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Figure  2.4.1 1.  Temperature  and  water  concentration  data  obtained  with  WMS-2/ showing  a  portion  of 
the  start-up  transient  after  combustor  ignition.  1392/1469  nm  line  pair.  <)>  =  1.20. 

Figure  2.4.12  shows  measurements  of  mean  temperature  as  a  function  of  fuel/air 
equivalence  ratio  using  the  1392/1469  nm  line  pair.  The  error  bars  represent  a  single 
standard  deviation  of  the  data  at  2  kHz.  Also  shown  is  the  temperature  calculated  by  a  1- 
D  model  developed  by  the  researchers  at  Wright-Patterson  AFB.  The  model  solves  the  1- 
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Figure  2.4.12.  Average  measured  temperature  and  1-D  simulated  temperature  (6.5  inches  downstream 
of  measurement  location)  versus  combustor  fuel/air  equivalence  ratio.  Error  bars  represent  standard 
deviation  at  2  kHz  sample  rate.  1 392/1469  nm  line  pair. 


D  conservation  equations  incorporating  reactant  mass  flow  rate,  load  cell  force,  heat  loss, 
and  pressure  data  from  each  run.  Some  deviation  between  the  simulated  and  measured 
temperatures  is  expected  because  the  model  predicts  the  temperature  at  the  end  of  the 
combustor,  which  is  16.5  cm  downstream  from  the  measurement  location.  Two 
competing  phenomena  occur  in  the  scramjet  between  the  measurement  location  and  the 
end  of  the  combustor.  The  duct  diverges,  causing  pressure  drop  and  a  corresponding 
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temperature  drop,  and  heat  release  still  may  be  occurring  as  the  gases  pass  the 
measurement  location,  which  will  cause  an  increase  in  temperature.  Another  factor  to 
consider  is  that  the  beam  path  for  the  laser  measurements  may  be  positioned  at  an 
elevation  that  is  not  representative  of  the  1-D  temperature.  The  key  aspect  of  Fig.  2.4.12 
is  that  the  trend  of  the  measured  temperature  matches  that  of  the  simulated  temperature. 
Measurements  of  this  nature  are  helpful  to  determine  optimal  scramjet  run  parameters 
such  as  fuel/air  ratio  or  to  characterize  the  efficiency  of  various  flame  holder  designs  or 
engine  components. 

Current  Status 

A  second-generation  scanned  WMS-2/  diode  laser  absorption  sensor  has  been 
designed,  built,  and  tested  in  the  laboratory  for  further  measurement  campaigns  in  Test 
Cell  22  at  the  AFRL  Propulsion  Directorate,  Wright-Patterson,  AFB,  OH.  This  system, 
termed  WMS-2///  incorporates  a  new  method  of  signal  normalization  that  accounts  for 
beam  steering,  window  fouling,  and  other  noise  sources.  In  addition,  this  normalization 
enables  quantitative  measurements  without  in  situ  calibration.  The  first  measurements 
using  this  new  sensor  system  are  planned  for  December  2006. 
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2.4.3.  Mid-IR  measurements  of  fuels  using  diode  lasers 

In  our  recent  AFOSR-sponsored  work  we  successfully  demonstrated  quantitative 
ethylene  and  propane  sensing  using  commercially  available  near-infrared  (NIR)  tunable 
diode  laser  sources  [2.4.6].  However,  hydrocarbon  absorption  at  these  NIR  wavelengths 
uses  transitions  in  the  vibrational  overtone  and  combination  bands  where  the  absorption 
cross  sections  are  relatively  small  and  that  become  even  smaller  and  more  problematic 
for  the  larger  and  more  complex  hydrocarbon  species  found  in  jet  fuels.  The  recent 
extension  of  solid-state  tunable  laser  technology  into  the  mid-infrared  (MIR)  now 
provides  the  opportunity  for  absorption  measurements  of  hydrocarbons  using  transitions 
in  the  fundamental  C-H  stretching  vibrational  mode  where  the  absorption  cross  section  is 
much  larger.  Using  this  much  stronger  absorption  strength,  new  and  much  more  sensitive 
diagnostic  techniques  can  be  anticipated.  However,  these  techniques  will  differ  from 
those  developed  previously  (in  the  NIR)  for  small  hydrocarbons  with  discrete  spectra 
owing  to  the  broader  spectral  features  of  relevant  large  hydrocarbon  molecules.  Recently 
we  have  begun  to  investigate  a  novel  dual-wavelength  differential  absorption  (DWDA) 
fuel  sensor  that  we  believe  holds  high  potential.  Our  first  proof-of-concept  experiments 
on  n-dodecane  are  described  below. 

Background 

The  enabling  technology  for  this  new  fuel  sensing  strategy  is  the  recent  commercial 
development  of  all-solid-state,  wavelength-tunable  mid-IR  laser  sources  using  difference- 
frequency-generation  (DFG)  and  our  acquisition  of  a  prototype  device  using  support  from 
our  AFOSR  Defense  University  Instrumentation  Program  (DURIP)  grant.  The  ability  to 
tune  or  select  the  specific  sensor  wavelength  is  very  important  for  mid-IR  absorption 
sensing.  Tunability  allows  selection  of  a  wavelength  that  avoids  potential  interference 
from  structured  absorption  of  water  vapor  combustion  products.  Wavelength  choice  also 
allows  the  absorption  strength  to  be  selected  from  the  broad  unstructured  hydrocarbon 
vibrational  absorption  features  for  maximum  dynamic  range  of  the  absorption  signal  for 
the  expected  range  of  fuel  concentration  and  path  length. 

The  spectroscopy  of  complex  hydrocarbon  fuel  species  is  significantly  different  than 
observed  for  small  polyatomic  compounds  such  as  H2O,  CH4,  C2H4,  or  C2H2,  which  all 
have  well-resolved  rotational-state-specific  structure.  The  rotational  structure  in  complex 
hydrocarbon  fuel  molecules  is  spaced  very  closely  and  blended  by  pressure  and  Doppler 
broadening  into  unresolved  absorption  features.  The  broad  structure  of  the  absorption 
spectrum  motivates  our  work  to  develop  a  differential  absorption  strategy  where 
measurements  are  made  at  two  colors  in  one  of  these  broad  absorption  features. 
Differential  absorption  measurements  will  allow  quantitative  measurement  of  the  vapor 
concentration  even  with  attenuation  of  the  laser  beam  from  window  fouling,  scattering 
losses,  or  beam  steering. 

Wavelength-tunable  mid-IR  laser  source  using  DFG 

The  fiber  optics  and  lasers  developed  for  telecommunications  in  the  NIR  have  the 
mechanical  robustness  needed  for  routine  sensing  in  aeroengine  research,  and  our  current 
AFOSR-sponsored  research  has  exploited  these  components  for  successful  measurements 
of  time-resolved  gas  temperature  in  scramjet  and  gas  turbine  combustors  [2.4.1]. 
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Recently  non-linear  conversion  techniques  have  extended  lasers  based  on  these  robust 
optical  components  into  the  mid-IR.  Wavelength-tunable  mid-IR  lasers  using  difference- 
frequency-generation  in  periodically  poled  lithium  niobate  (PPLN)  have  been  used  for 
the  past  decade  for  sensitive  detection  of  trace  species  in  the  atmosphere,  but  virtually  all 
of  this  work  has  required  use  of  large  pulsed  lasers.  The  availability  of  continuous  (cw) 
sources  has  remained  an  elusive,  unfulfilled  development  goal. 

Recently,  however,  a  cw  source  using  tunable  diode  lasers  and  fiber  amplifiers  was 
demonstrated,  allowing  reduced  power  consumption  for  balloon-borne  measurements. 
This  laser  design  proved  quite  stable  and  robust,  and  a  similar  architecture  now  is  being 
commercialized  by  a  San  Francisco  start-up  company  (NovaWave).  The  laser  design, 
illustrated  in  Fig.  2.4.13,  fiber-amplifies  a  NIR  wavelength-tunable  diode  laser  (seed)  and 
combines  this  light  in  a  polarization-maintaining  fiber  with  a  powerful  single-frequency 
diode  laser.  These  two  near-IR  laser  wavelengths  are  mixed  in  a  periodically  poled 
lithium  niobate  (PPLN)  crystal  to  produce  laser  light  at  the  difference  frequency  (here 
near  3.4  pm  in  the  mid-IR).  The  mid-IR  light  can  be  tuned  in  wavelength  by  changing 
the  wavelength  of  the  near-IR  seed  laser,  and  when  this  change  is  done  by  rapidly  (thus 
far  we  have  demonstrated  up  to  100  kHz)  alternating  between  two  separate  near-IR  seed 
lasers,  as  illustrated  in  Fig.  2.4.13,  the  output  consists  of  a  beam  that  alternates  between 
two-colors  of  mid-IR  light  along  a  single  optical  path.  The  current  laser,  recently 
procured  on  our  AFOSR  DURIP  grant  and  based  on  a  wavelength-dithering  concept  we 
suggested  to  the  vendor,  enables  us  to  separate  the  two  mid-IR  colors  by  ~10  cm'1,  which 
is  well-suited  for  the  initial  feasibility  measurements  described  below  on  evaporating  n- 
dodecane  vapor.  Custom  PPLN  crystals  would  enable  this  color  separation  to  be 
enhanced  significantly,  enabling  further  advances  in  sensor  performance  for  cases  where 
increased  spectral  separation  would  yield  better  vapor  sensing  capability  in  the  presence 
of  liquid  fuel  droplets. 
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Figure  2.4.13.  Architecture  of  all  solid-state  two-color  DFG  mid-IR  light  source. 
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Dual-wavelength  differential  absorption  to  suppress  aerosol  interference 

We  demonstrated  the  dual-wavelength  differential-absorption  (DWDA)  sensor 
concept  using  this  novel  two-color  mid-infrared  (mid-IR)  laser.  The  total  Beer’s  law 
extinction  at  each  frequency  v  (color)  is  the  sum  of  vapor  absorption  and  droplet 
scattering: 

—  exp[-  ccvap0r  (>')—  T  droplets  (v)] 
v 
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Thus  if  we  measure  the  extinction  at  two  laser  frequencies,  V|  and  V2,  solving  for  the 
vapor  absorption  reduces  to  the  solution  of  two  equations  with  two  unknowns.  This 
DWDA  sensing  concept  is  able  to  extract  the  vapor  absorption  signal  from  strong 
interference  attenuation  due  to  scattering  from  the  aerosol  droplets  during  evaporation  of 
shock-heated  n-dodecane  aerosol.  First  we  examine  the  absorption  spectrum  to  select  the 
laser  frequencies.  The  left  panel  of  Fig.  2.4.14  shows  the  mid-IR  absorption  spectrum  of 
n-dodecane  in  the  region  of  the  C-H  stretching  vibration.  We  select  two  laser  frequencies, 
V|=2916cm'1  and  V2=2926cm'',  to  maximize  the  difference  in  the  absorption  coefficient 
between  the  two  laser  colors,  as  illustrated  on  the  plot.  Then  we  must  understand  the 
dependence  of  the  absorption  spectrum  on  temperature  (and  pressure).  The  measured 
pressure  dependence  is  small  and  can  be  neglected  for  n-dodecane  for  P<10atm; 
however,  there  is  a  significant  temperature  dependence,  as  illustrated  in  the  right  panel  of 
Fig.  2.4.14  where  the  differential  absorption  coefficient  is  shown  versus  T. 
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Figure  2.4.14.  Left  panel:  absorption  coefficient  of  n-dodecane  vapor  dilute  in  1  atmosphere  of 
nitrogen  measured  at  675K  via  FTIR;  Right  panel:  differential  absorption  coefficient  versus  gas 
temperature  for  n-dodecane  vapor  at  2916  and  2926  cm'1. 

An  aerosol  of  n-dodecane  is  prepared  via  an  ultrasonic  nebulizer  and  mixed  with 
argon  in  Stanford  University’s  aerosol  shock  tube  facility,  as  illustrated  in  Fig.  2.4.15. 
DWDA  of  the  two-color  mid-IR  laser  extinction  is  used  to  extract  the  vapor 
concentration  of  n-dodecane,  and  extinction  of  the  NIR  laser  is  used  to  monitor  the  vapor 
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Figure  2.4.15.  Stanford’s  aerosol  shock  tube  with  the  mid-IR  DWDA  diagnostic  plus  NIR  extinction  to 
monitor  the  droplet  loading. 
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loading.  The  two-color  mid-IR  laser  switches  between  V|=2916cm'1  and  V2=2926cm''  at 
100  kHz,  thus  enabling  differential  absorption  measurements  with  10  ps  time  resolution. 
The  measurement  is  shown  in  Fig.  2.4.16  for  a  pre-shock  pressure  of  0.22  atm  and  pre¬ 
shock  temperature  of  300K.  Before  the  shock-heating,  the  aerosol  scattering  reduces  the 
laser  transmission  by  48%,  as  seen  by  the  NIR  transmission.  Even  with  this  significant 
aerosol  scattering  interference,  the  DWDA  measurement  recovers  the  small  300K  vapor 
pressure  (Pvapor=0.16  Torr)  of  n-dodecane.  When  the  incident  shock  arrives,  the  aerosol- 
laden  gas  is  compressed  to  a  post-shock  pressure  of  0.78  atm  and  post-shock  temperature 
of  436K.  The  aerosol  loading  increases  because  of  shock  compression,  and  the 
transmission  of  the  NIR  laser  light  drops  below  20%  (i.e.,  the  interference  extinction  is 
greater  than  80%!),  as  seen  in  Fig.  2.4.16.  However,  at  the  post-shock  temperature  of 
436K,  the  aerosol  rapidly  evaporates,  and  the  time-resolved  increase  in  the  n-dodecane 
vapor  is  monitored  successfully  by  DWDA,  as  shown  in  the  figure.  When  the  droplets 
are  evaporated  completely  the  NIR  transmission  recovers  to  100%,  and  the  DWDA 
measures  the  vapor  initially  contained  in  the  droplets  plus  the  initial  300K  vapor  pressure. 


Figure  2.4.16.  Shock-heated  n-dodecane  aerosol  monitored  via  mid-IR  DWDA  for  vapor  mole 
fraction  of  do-decane  and  NIR  transmission  to  monitor  aerosol  loading:  Test  gas  at  300K, 
P=0.22  atm,  Xv«por=0.09%  with  post-shock  T=436K,  P=0.78atm  and  final  X«poc=0.54%. 


Current  Status 

The  measurements  shown  in  Fig.  2.4.16  formed  the  cornerstone  of  a  presentation  at 
the  31st  International  Combustion  Symposium  last  summer  [2.4.7].  We  believe  this  new 
concept  for  diode  laser-based  mid-infrared  sensing  offers  significant  potential  for  sensing 
hydrocarbon  fuels  for  fundamental  chemistry  experiments,  laboratory  combustion 
experiments,  and  practical  aeroengine  ground  test. 
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2.5  Photophysics  studies  in  support  of  tracers  for  fuel  imaging 

Background 

PLIF  imaging  of  tracer  molecules  has  proven  to  be  a  powerful  diagnostic  for 
investigating  the  behavior  of  fuels,  for  example  in  mixing  and  combustion.  PLIF 
imaging  of  acetone,  pioneered  at  Stanford  with  AFOSR  support,  is  widespread  owing  to 
its  ease  of  use  and,  importantly,  to  its  simple  interpretation  for  isobaric  and  isothermal 
flows.  Similarly,  3-pentanone,  a  slightly  larger  member  of  the  ketone  family,  has  seen 
growing  use,  particularly  in  internal  combustion  engines.  Past  work  at  Stanford  has 
sought  build  a  fundamental  understanding  of  the  photophysics  of  these  tracers  to  enable 
their  quantitative  use  in  flows  with  varying  pressure,  temperature,  and  composition. 
Accurate  databases  for  the  absorption  and  fluorescence  behavior  of  these  two  ketones 
have  thus  been  acquired  and  physically  realistic  photophysical  models  assembled  and  fit 
to  the  available  experimental  data.  These  models,  in  turn,  have  been  used  to  identify  new 
diagnostic  strategies  for  simultaneously  measuring  the  tracer  concentration  and 
temperature,  as  well  as  to  enable  quantitative  application  of  PLIF  imaging  over  a  range  of 
pressures  and  temperatures.  This  body  of  work  on  tracer  photophysics  and  modeling  has 
had  considerable  impact  internationally  on  combustion  diagnostics. 

Current  Status  of  Toluene  Photophvsics 

During  the  past  AFOSR-grant  period,  we  have  begun  to  extend  this  work  on  the 
photophysics  of  fuel  tracers  to  include  toluene,  an  aromatic  compound  that  is  present  in 
jet  fuel  and  is  commonly  selected  as  the  aromatic  component  in  surrogate  mixtures  of 
hydrocarbons  used  to  simulate  behavior  of  jet  fuels.  While  we  are  still  in  a  relatively 
early  phase  of  these  studies  with  toluene,  we  already  have  acquired  a  database  for  the 
absorption  spectrum  as  a  function  of  temperature  and  for  the  fluorescence  yield  at  two 
selected  excitation  wavelengths  (248  nm  and  266  nm)  but  essentially  limited  to  nitrogen 
as  a  bath  gas.  This  work  has  documented  a  very  strong  temperature  dependence  of  the 
toluene  LIF  signal,  shown  in  Fig.  2.5.1,  which  suggests  exciting  potential  for  a  new, 
highly  sensitive  PLIF  temperature  diagnostic  [2.5.1].  There  is  also  a  very  strong 
dependence  of  the  toluene  fluorescence  on  oxygen  concentration,  which  leaves  open  the 
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Figure  2.5.1.  Relative  toluene  LIF  signal  from  248-nm  and  266-nm  excitation  as  a  function  of 
temperature  at  1  bar  total  pressure  (diluted  with  nitrogen).  [2.5.1] 
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hope  of  simultaneously  imaging  temperature  and  fuel-air  ratio  [2.5.2],  As  an  illustrative 
example  of  the  promise  of  a  toluene  PLIF  diagnostic.  Fig.  2.5.2  provides  a  single-shot,  2- 
D  temperature  measurement  image  acquired  for  a  heated  cylinder  in  a  uniform  nitrogen 
flow  seeded  with  toluene.  The  high  precision  of  the  measurement  enables  a  resolution  of 
temperature  differences  of  about  5  K. 


Figure  2.5.2.  Single-shot  PLIF  temperature  image  of  uniform  nitrogen  flow  around  a  heated 
cylinder  using  a  toluene  tracer.  The  laser  sheet  is  incident  from  the  left  side  of  the  image. 


Research  into  High-Temperature  and  -Pressure  Toluene  Photophvsics 

To  characterize  the  PLIF  signal,  we  utilize  the  fluorescence  equation: 

Sf  oc  E  n-  a(A,  T)  •  <f>(A,  T,  p ) 

where  E  is  laser  energy,  n  is  the  tracer  number  density,  a  is  the  absorption  cross  section, 
and  <j>  is  the  fluorescence  quantum  yield.  Any  quantitative  PLIF  diagnostic  clearly 
depends  on  knowledge  of  the  absorption  cross  section  and  the  fluorescence  quantum 
yield.  Measurements  of  the  variations  in  the  absorption  cross  section  and  the 
fluorescence  quantum  yield  with  respect  to  pressure,  temperature,  and  composition  are 
thus  of  special  importance  to  the  development  of  a  quantitative  and  optimized  toluene 
PLIF  diagnostic,  and  very  little  data  of  this  type  currently  exist  for  toluene. 

As  stated  previously,  our  preliminary  measurements  of  the  fluorescence  quantum 
yield  in  toluene  and  its  variation  with  temperature  are  conducted  in  nitrogen  bath  gas  at  a 
pressure  of  1  atm.  The  experimental  setup  used  to  make  these  measurements  is  shown  in 
Fig.  2.5.3,  and  recently  we  have  designed  modifications  to  this  test  cell  to  accommodate  a 
much  wider  range  of  pressures  and  temperatures. 

The  tracer  (e.g.  toluene)  is  introduced  into  the  carrier  gas  using  a  saturating 
percolator.  The  seeded  gas  subsequently  flows  into  an  optically  accessible  cell,  where  it 
is  excited  by  the  appropriate  laser.  The  fluorescence  signal  is  measured  via  either  a  PMT 
or  a  CCD  camera  mounted  on  a  spectrograph.  By  using  a  furnace  with  the  cell,  the  effect 
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Figure  2..5J.  Experimental  setup  for  measuring  absorption  cross  section  and  fluorescence  quantum 
yield  of  tracers  as  a  function  of  temperature. 

of  increasing  temperature  on  the  fluorescence  signal  can  be  measured.  By  increasing  the 
internal  pressure  of  the  cell,  the  effect  of  pressure  of  fluorescence  signal  can  be  measured. 
This  measurement  technique,  including  a  novel  Rayleigh  scattering  approach  for  absolute 
calibration  of  fluorescence  yield,  is  developed  in  work  on  ketone  fluorescence  [2.5.3]  and 
has  worked  well  in  our  initial  work  with  toluene. 

Semi-Empirical  Fluorescence  Model  for  Ketone  and  Aromatic  Tracers 

A  validated,  physics-based  model  for  the  fluorescence  quantum  yield  is  critical  for 
the  development  of  successful  diagnostic  strategies.  The  fluorescence  signal  can  be 
characterized  using  the  fluorescence  equation: 

Sf  oc  E ■  n ■  a(A,T)  <f>(A,T,p) 

The  critical  information  needed  is  embodied  in  the  absorption  cross  section  and  the 
fluorescence  quantum  yield.  A  model  of  fluorescence  quantum  yield  not  only  will 
contribute  to  the  understanding  of  the  processes  underlying  PL1F,  but  also  will  allow  for 
the  prediction  of  fluorescence  behavior  and  optimization  of  diagnostic  strategies  for 
various  conditions.  Our  past  work  to  model  the  fluorescence  quantum  yield  for  acetone 
and  3-pentanone  has  contributed  to  the  expanding  use  of  ketone  fluorescence  in  the  PLIF 
diagnostic  arena  [2.5.4]. 

Fig.  2.5.4  diagrams  our  current  conceptualization  of  the  fluorescence  process 
developed  for  ketones  and  extended  for  toluene.  The  laser  photon  excites  the  molecule 
from  the  So  state  to  the  Si  state.  In  the  Si  state,  in  addition  to  fluorescence,  the  molecule 
could  undergo  several  processes.  Intramolecular  vibrational  redistribution  (IVR) 
redistributes  the  vibrational  energy  among  all  the  levels  of  S|  without  changing  total 
vibrational  energy.  IVR  also  affects  intersystem  crossing  (ISC)  and  internal  conversion 
(IC)  rates,  as  these  rates  depend  on  the  vibrational  energy  in  the  S|  state.  In  ISC,  the 
molecule  goes  from  the  excited  state  S|  to  the  triplet  state  T|.  In  IC,  the  molecule  goes 
from  the  Si  state  back  to  the  So  state  via  the  upper  vibrational  levels  of  the  So  state.  In 
addition,  there  is  potential  for  collisions  with  oxygen  molecules  to  promote  further 
transfer  of  the  toluene  molecule  to  the  So  and  T|  states  through  formation  of  short-lived 
molecular  complexes.  Besides  fluorescence,  none  of  these  processes  emit  photons,  so  the 
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sole  source  of  signal  is  fluorescence.  The  variation  of  these  competing  processes  with 
pressure,  temperature,  and  bath  gas  composition  holds  the  key  to  allowing  use  of  the 
fluorescence  signal  as  a  quantitative  gasdynamic  diagnostic. 


1C  ISC 


Figure  2. .5.4.  Fluorescence  energy  diagram  for  toluene. 
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Pressure  Combustion  Systems,”  Ph.D.  Thesis,  Department  of  Mechanical 
Engineering,  Stanford,  CA,  March,  2006;  Thermosciences  Division  Report  1 65, 
http://navier.stanford.edu/thermosciences/TSD-165.pdf. 

7.  Xiang  Liu,  “Line-of-Sight  Absorption  of  H2O  Vapor:  Gas  Temperature  Sensing  in 
Uniform  and  Nonuniform  Flows,”  Ph.D.  Thesis,  Department  of  Mechanical 
Engineering,  Stanford,  CA,  June,  2006;  Thermosciences  Division  Report  1 68, 
http://navier.stanford.edu/thermosciences/TSD- 1 68.pdf. 
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3.3  AFOSR  SPONSORED  PRESENTA  TIONS  2004-2006 

MEETING  PRESENTATIONS  2004 


I .  Two  presentations  at  the  AIAA  Aerospace  Science  Meeting  (papers  listed  above) 

3.  One  presentation  at  the  AIAA  AMT  and  Ground  Test  Meeting  (paper  listed  above) 

4.  Three  oral  presentations  at  the  30th  International  Symposium  on  Combustion  in 
Chicago  (reviewed  papers  listed  as  published  in  2005) 

7.  W.  Koban,  J.D.  Koch,  R.  K.  Hanson,  and  C.  Schulz,  “FARL1F:  Myths  and  Reality,”  poster 
at  the  3(fh  International  Symposium  on  Combustion  in  Chicago,  July  2004. 

8.  W.G.  Bessler,  H.  Kronemayer,  C.  Schulz,  T.  Lee,  J.B.  Jeffries,  R.K.  Hanson, 
“Quantitative  multi-line  NO-LIF  temperature  imaging  in  flames  over  a  wide 
pressure  range,''’  Laser  Applications  to  Chemical,  Environmental,  and  Atmospheric 
Analysis,  Optical  Society  of  America,  Annapolis,  Feb.  2004. 

9.  E.  Barbour,  M.A.  Oelschlaeger,  D.W.  Mattison,  D.F.  Davidson,  C.  Schulz,  J.B. 
Jeffries,  R.K.  Hanson,  “UV  absorption  of  CO2  for  temperature  diagnostics,”  Laser 
Applications  to  Chemical,  Environmental,  and  Atmospheric  Analysis,  Optical 
Society  of  America,  Annapolis,  Feb.  2004. 

1 0.  J.B.  Jeffries  (invited),  “Tunable  diode  laser  sensors  for  practical  combustion 
applications,”  Conference  on  Lasers  and  Electro-Optics,  Optical  Society  of 
America,  San  Francisco,  CA,  May,  2004. 

MEETING  PRESENTATIONS  2005 

I I .  One  presentation  at  the  AIAA  Aerospace  Sciences  Meeting  (paper  listed  above) 

12.  Two  presentations  at  the  AIAA  JPC  Conference  (papers  listed  above) 

13.  D.W.  Mattison,  J.T.C.  Liu,  G.B.  Rieker,  J.B.  Jeffries,  R.K.  Hanson,  “Diode  laser 
sensors  for  internal  combustion  and  aeropropulsion  engine  applications,” 
Thermosciences  Division  Industrial  Affiliates,  Stanford  University,  February,  2005. 

14.  T.  Lee,  W.G.  Bessler,  H.  Kronemayer,  C.Schulz,  J.B.  Jeffries,  R.K.  Hanson, 
“Quantitative  temperature  measurements  in  high-pressure  flames  with  multi-line 
nitric  oxide  (NO)-LIF  thermometry,”  Joint  Meeting  of  the  US  Sections  of  the 
Combustion  Institute,  Philadelphia,  PA,  March,  2005. 

15.  R.K.  Hanson,  J.B.  Jeffries,  M.G.  Allen,  “Tunable  diode  laser  absorption  sensor 
applications  to  aeropropulsion  testing,”  NA TO  RSM-01 7  Applied  Vehicle 
Technology  Panel  Specialists  ’  Meeting  on  Recent  Developments  in  Non-Intrusive 
Measurement  Technology,  Budapest,  Hungary,  April,  2005;  invited  plenary  lecture. 

16.  M.  Gruber,  D.  Eklund,  T.  Mathur,  J.  Jeffries,  "In-stream  measurements  at  the  exit  of 
dual-mode  hydrocarbon-fueled  scramjet  combustors,"  4Cfh  Combustion 
Subcommittee,  2&h  Airbreathing  Propulsion  Subcommittee ,  22nd  Propulsion 
Systems  Hazards  Subcommittee,  and  4,h  Modeling  and  Simulation  Subcommittee 
Joint  Meeting,  JANNAF,  June  2005. 
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1 7.  R.K.  Hanson  (invited)  and  J.B.  Jeffries,  “Advances  in  diode  laser  sensors  for 
combustion  and  propulsion,”  Western  States  Section  Meeting  of  the  Combustion 
Institute ,  paper  05-01 F  Stanford,  CA,  October,  2005. 

MEETING  PRESENTATIONS  2006 

1 8.  One  presentation  at  the  AIAA  Aerospace  Sciences  Meeting  (paper  listed  above) 

19.  One  presentation  at  the  AIAA  AMT  and  Ground  Test  Meeting  (paper  listed  above) 

20.  Two  presentations  at  the  AIAA  JPC  Conference  (papers  listed  above) 

22.  Three  presentations  at  the  31st  International  Symposium  on  Combustion  in 
Heidelberg  Germany  (reviewed  papers  listed  as  published  in  2007) 

26.  Ronald  K.  Hanson  (invited)  and  Jay  B.  Jeffries,  “Advances  in  diode  laser  absorption 
sensors  for  combustion  and  propulsion,”  Laser  Applications  to  Chemical, 
Environmental,  and  Atmospheric  Analysis,  Optical  Society  of  America,  Incline 
Village,  CA,  February,  2006. 

27.  R.K.  Hanson  and  J.B.  Jeffries,  “Diode  sensors  for  combustion  and  propulsion,”  41st 
Combustion  Subcommittee,  2&h  Airbreathing  Propulsion  Subcommittee ,  23rd 
Propulsion  Systems  Hazards  Subcommittee,  JANNAF,  San  Diego,  CA,  December, 
2006. 

INVITED  LECTURES  2004-2006 

28.  R.K.  Hanson,  “Advances  in  planar  laser-induced  fluorescence  imaging  and 
applications  to  combustion  and  propulsion,  JAXA,  Japan  Space  Agency,  March, 
2004. 

29.  R.  K.  Hanson,  “Tunable  diode  lasers  for  sensing  and  control  with  applications  in 
combustion  and  propulsion,  Tokyo  Institute  of  Technology,  March  2004. 

30.  J.B.  Jeffries,  “Laser  diagnostics  for  combustion  and  propulsion.  University  of 
Duisburg-Essen,  Duisburg,  Germany,  April  2004 

31.  R.K. Hanson  and  J.B.Jeffries,  “Advances  in  laser-based  sensors  for  combustion 
research,”  Sixth  Symposium  on  Smart  Control  of  Turbulence,  Tokyo,  Japan,  March, 
2005;  invited  plenary  lecture. 

32.  R.K.  Hanson,  “Diode  laser  sensors  for  combustion  and  propulsion  applications,” 
Tokyo  Institute  of  Technology,  March,  2005. 

33.  M.  A.  Oehlschlaeger,  “Shock  tube  studies  of  thermal  decomposition  reactions  using 
ultraviolet  laser  absorption  spectroscopy,”  Combustion  Research  Facility,  Sandia 
National  Laboratory,  Livermore,  CA,  May,  2005. 

34.  J.  Koch,  “Fuel  tracer  photophysics  for  quantitative  planar  laser-induced 
fluorescence,”  Marquette  University,  April,  2005. 

35.  J.B.  Jeffries,  “UV  absorption  of  CO2:  opportunities  for  new  combustion 
diagnostics,”  Gordon  Research  Conference  on  Laser  Diagnostics  for  Combustion, 
August,  2005. 
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36.  R.K.  Hanson,  “Optical  Diagnostics  for  Hypersonic  Flows,”  University  of 
Minnesota/ AFOSR  Workshop  on  Hypersonics,  September,  2005. 

37.  R.K.  Hanson,  “TDL  Sensing  for  Augmentor  Applications,”  Briefing  at  WPAFB 
November,  2005. 

38.  Ronald  K.  Hanson,  “Laser-based  Sensing  and  Control  of  Combustion 
Applications,”  Workshop  on  Control  and  Enhancement  of  Combustion  of  High 
Velocity  Vitiated  Flows  by  Pulsed  Discharges,  WPAFB,  November,  2005. 

39.  Jay  B.  Jeffries  and  Ronald  K.  Hanson,  “Advances  in  Diode  Laser  Sensors  for 
Propulsion  and  Combustion,  AFRL,  WPAFB,  March  2006. 

40.  Ronald  K.  Hanson,  “Advances  in  Tunable  Diode  Laser  (TDL)  Sensing  for 
Combustion  and  Propulsion,”  Keynote  Address,  Georgia  Tech  Symposium  in  Honor 
of  Professor  Ben  Zinn,  May,  2006. 

41.  Jay  B.  Jeffries  and  Ronald  K.  Hanson,  “Optical  Diagnostics  for  Practical 
Applications,”  Instrumentation  Roundtable  on  Ballistic  Fire  Initiation,  Joint 
Aircraft  Survivability  Program  Review,  Nellis  AFB,  September,  2006. 

42.  Ronald  K.  Hanson,  “Applications  of  Laser  Diagnostics  to  Combustion  and 
Propulsion,”  Lund  University  Laser  Center,  Lund,  University,  September,  2006. 

43.  Ronald  K.  Hanson,  “Advances  in  Diode  Laser  Sensors  for  Combustion  and 
Propulsion,  Mechanical  Engineering  Seminar,  Purdue  University,  September,  2006. 
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4.0  Personnel 


Individual  researchers  partially  or  fully  supported  by  the  program  during  the 
reporting  period  are  listed  below.  All  the  work  has  been  carried  out  in  the  High 
Temperature  Gasdynamics  Laboratory,  in  the  Department  of  Mechanical  Engineering, 
under  the  supervision  of  Professor  R.  K.  Hanson. 

Research  Staff 

Dr.  Jay  B.  Jeffries,  Senior  Research  Associate 

Current  Graduate  Research  Assistants 

Brian  Chung 
Andrew  Fahrland 
Hejie  Li 
Gregory  Rieker 
David  Rothamer 
John  Yoo 

Ph.D  Degrees  Awarded  2004-2006  with  Partial  AFOSR  Support 

Jonathan  Liu,  2004,  currently  NIH  fellow,  Stanford  University 

Jon  Koch,  2005,  currently  Assistant  Professor,  Marquette  University 

M.  A.  Oehlschlaeger,  currently  Assistant  Professor,  Rensselaer  Polytechnic  Inst. 

Lin  Ma,  2005,  currently  Assistant  Professor,  Clemson  University 

Dan  Mattison,  2006,  currently  with  Exponent,  Inc. 

Tonghun  Lee,  2006,  currently  Assistant  Professor,  Michigan  State  University 
Xiang  Liu,  2006,  currently  with  Spectral  Sensors,  Inc. 

Visitors  Supporting  the  AFOSR  Effort 

Dr.  Christof  Schulz,  2004,  Visiting  Professor  from  the  University  of  Heidelberg 
Wieland  Koban,  2004,  Visiting  Scientist  from  the  University  of  Heidelberg 

Wolfgang  Bessler,  2004,  Visiting  Scientist  from  the  University  of  Heidelberg 
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5.0  Significant  Interactions 


In  addition  to  the  interactions  associated  with  the  presentations  and  publications 
listed  in  Section  3  we  have  had  numerous  visitors  to  our  laboratory  during  the  contract 
period.  Foreign  visitors  have  come  from  the  United  Kingdom,  Germany  and  Japan; 
industrial  and  national  laboratory  visitors  have  included  representatives  from: 

Industry: 

Aerometrics 

Boeing 

Exxon 

NEL  and  NTT 

GE  Research  Laboratories 

Metrolaser 

Physical  Sciences,  Inc. 

Pratt  and  Whitney 
Rolls  Royce 
TRW 

United  Technologies 
Zolo  Technologies 
Government 

Air  Force  Research  Laboratory  (WPAFB) 

Army  Research  Laboratories  (ARL) 

NASA  Ames 
NASA  Glen 

Sandia  Livermore  National  Laboratory. 

5.7  Interactions  and  Collaboration  with  AFRL 

Stanford  University  hosted  a  workshop  “Understanding  Plasma  Ignition”  on  January 
9,  2004.  Attendees  included  AFRL  and  AFOSR  personnel;  Skip  Williams, 
AFRL/Geophysics;  Campbell  Carter,  AFRL/Propulsion;  Julian  Tiskoff,  AFOSR;  and 
Michael  Berman,  AFOSR.  AFOSR  contractors  were  also  participants  including: 
Svetlana  Starikovskaia,  MIPT;  Andrei  Starikovskii,  MIPT;  Vincent  McCoy,  Cal  Tech, 
Mark  Kushner,  U.  Illinois;  Robert  Continetti,  UCSD;  Martin  Gundersen,  USC;  Jay 
Jeffries  and  Ronald  Hanson,  Stanford  University.  In  addition,  one  participant,  Jose 
Sinibaldi,  came  from  the  Naval  Postgraduate  School  in  Monterey. 

We  collaborated  with  Dale  Shouse  of  AFRL  at  WPAFB  during  2003  on  the  writing  of 
a  paper  describing  the  diode  laser  sensor  measurements  performed  in  the  gas  turbine 
sector  test  rig  at  AFRL  in  November  of  2002.  Jonathan  Liu  of  Stanford  University 
presented  this  paper  at  the  AIAA  Joint  Propulsion  Conference  in  Huntsville  Alabama  in 
July  of  2003  and  won  AIAA’s  Ground  Test  best  paper  award  in  June  of  2004 

A  measurement  campaign  was  conducted  in  September  2004  to  perform  diode  laser 
sensing  measurements  in  the  scramjet  combustor  test  facility  at  AFRL  at  Wright 
Patterson  AFB,  OH.  This  was  the  second  set  of  measurements  conducted  in 
collaboration  with  Dr.  Mark  Gruber  and  Dr.  Campbell  Carter  of  AFRL.  The  improved 
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optical  engineering  used  to  match  this  sensor  with  the  optical  access  in  the  test  facility 
enabled  major  improvement  in  signal  to  noise.  The  results  have  been  presented  in  at  the 
AIAA  Joint  Propulsion  Conference  by  Stanford  student  Greg  Rieker  (Publications  and 
AIAA  papers  #21),  in  the  archival  journal  Applied  Optics  (Publication  #25),  and  included 
in  the  2005  JANNAF  presentation  by  Dr.  Mark  Gruber  of  AFRL  (Presentations  at 
Meetings  #16).  A  third  measurement  campaign  to  AFRL  to  investigate  this  measurement 
technology  for  H2O  concentration  and  gas  temperature  is  planned  for  December  2006. 
The  third  campaign  has  three  goals:  (1)  compare  the  direct  absorption  sensing  strategy 
used  previously  with  wavelength  modulation  spectroscopy  with  2f  detection,  (2) 
transition  from  experiments  to  investigate  sensor  performance  for  gas  temperature  to 
experiments  to  use  gas  temperature  to  study  the  performance  of  the  scramjet  combustor, 
and  (3)  investigate  the  feasibility  of  measurements  of  CO2  near  2pm. 

Dr.  Jeffries  of  Stanford  participated  in  a  workshop  in  May  2005  at  AFRL  organized 
by  Barry  Kiel  to  facilitate  diagnostic  measurements  on  augmented  jet  engines.  A  joint 
project  with  Pratt  and  Whitney  (Jeff  Lovett)  and  AEDC  has  begun  to  develop  a  prototype 
sensor  strategy  for  measurements  in  the  augmentor  of  high-performance  military  gas 
turbine  aeroengines. 

During  the  grant  period,  the  Hanson  group  had  significant  interaction  with  AFRL 
scientists  at  Wright-Patterson  AFB.  For  example,  in  addition  to  the  collaborative 
measurements  in  Test  Cell  22  at  WPAFB,  Professor  Hanson  participated  in  a  workshop  at 
in  November,  2005  on  plasma  stabilized  combustion.  Dr.  Jeffries  gave  a  public  talk  at 
AFRL  in  March  2006  on  diode  laser  sensors  for  propulsion  ground  test. 

During  this  period  of  AFOSR  support.  Professor  Hanson  has  given  invited 
presentations  to  several  industrial  laboratories,  universities,  and  government  groups  that 
highlight  our  progress  on  AFOSR-sponsored  diagnostics  research.  Members  of  our 
group  have  also  provided  technical  information  and  advice,  by  telephone  and  mail,  to 
several  external  researchers  interested  in  duplicating  or  extending  our  diagnostics 
concepts. 

5.2  Transitions 


Significant  technical  transitions  of  the  combustion  diagnostics  developed  under 
this  program  have  occurred  with  several  industrial  and  government  scientists  during  the 
grant  period.  This  work  has  primarily  focused  on  the  development  of  practical  sensors 
using  the  tunable  diode  laser  diagnostic  strategies  developed  over  the  past  decade  at 
Stanford  University.  These  significant  technology  transitions  include: 


Customer 

Result 

Application 

Dr.  Jeffery  Lovett, 

Pratt  and  Whitney, 

East  Hartford,  CT 
(860)557-0559 

Tunable  diode  laser 
diagnostics  for  gas 
temperature  measurements 

Diagnostics  of  instabilities  in 
augmented  aircraft  engines 
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Dr.  Michael  Aamio, 
Pratt  and  Whitney, 
Seattle  Aerosciences 
Center 

(425-450-3930) 

Tunable  diode  laser 
diagnostics 

Use  in  pulse  detonation 
engine  development  tests  at 
the  Naval  Warfare  Center  at 
China  Lake,  CA 

Dr.  Michael  Winter, 
Pratt  and  Whitney, 
East  Hartford  (860- 
565-9007) 

Tunable  diode  laser 
diagnostics 

Development  of  new  sensor 
for  air  mass  flux  for  in-flight 
aero-engine  testing;  proof-of- 
concept  tests  on  PW6000 
aeroengine  July  2005 

Dr.  Mark 
Woodmansee, 

General  Electric 

Global  Research 
Center  (518-387- 
4220) 

Tunable  diode  laser 
diagnostics 

Potential  new  sensor  for 
exhaust  pattern  factor  from 
fixed  gas  turbines  for  power 
generation 

Dr.  Anthony  Dean, 
General  Electric 

Global  Research 
Center  (518-387- 
6478) 

Tunable  NIR  and  MIR  laser 
diagnostics 

Stanford  demonstration  at  GE 
research  laboratories  in  pulse 
detonation  engine 
development  tests 

Dr.  Mark  Allen 

PSI  Corp 

Andover,  MA 
(978)689-0003 

Gas  sensors  based  on 
tunable  diode  lasers  using  2f 
techniques 

Temperature  sensors  for 
practical  industrial 
applications;  sensor  for  IC- 
engines  developed 

Dr.  Andrew  Sappey 
Zolo  Technologies 
Boulder,  CO 
(303)604-5800 

Wavelength-Multiplexed 
tunable  diode  laser  sensor 
strategies 

Temperature  sensors  for 
practical  industrial 
applications 

Dr.  Douglas  Baer, 
Los  Gatos  Research 
(650  965  7778) 

High  Pressure  Laser- 
Induced  Fluorescence  of 

NO 

Gas  turbine  combustor 
ground  test  measurements 
of  NO  spatial  distributions 

Dr.  Andrew 

Sappey,  Zolo 
Technology,  (303- 
604-5804 ) 

Tunable  diode  laser 
diagnostics 

Ground  test  of  augmented 
jet  engines  and  scramjet 
vehicle  health  monitoring 
systems 

6.0  New  discoveries,  inventions,  patent  disclosures 

No  patent  or  invention  disclosures. 


71 


7.0  Honors/ Awards 

Professor  Ronald  K.  Hanson  was  elected  Fellow  of  the  ASME  during  2004. 

The  AIAA  Best  Paper  Award  for  Ground  Testing  was  presented  at  the  40th  JPC  in 
Portland,  OR,  in  June,  2004  to  the  paper  by  J.T.C.  Liu,  J.B.  Jeffries,  R.K.  Hanson,  S. 
Creighton,  J.A.  Lovett,  D.T.  Shouse  entitled  “Diode  Laser  Absorption  Diagnostics  for 
Measurements  in  Practical  Combustion  Flow  Fields.”  This  paper,  which  had  been  given  at 
the  39th  AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference,  Huntsville,  Al,  July,  2003, 
included  authors  from  Stanford  University,  Pratt  and  Whitney,  and  the  Air  Force 
Research  Laboratory. 

Lifetime  awards  of  Professor  Ronald  K.  Hanson: 

1990  Fellow,  Optical  Society  of  America  (OSA) 

Woodard  Chair  of  Mechanical  Engineering,  Stanford  University 

1996  Aerodynamic  Measurement  Technology  award  from  the  AIAA  (first  recipient) 

1997  Fellow,  American  Institute  of  Aeronautics  and  Astronautics  (AIAA) 

2001  Elected  to  the  National  Academy  of  Engineering 

2002  Silver  Medal  of  the  Combustion  Institute 

2004  Fellow,  American  Society  of  Mechanical  Engineers 

2005  Propellants  and  Combustion  award  from  the  AIAA 
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